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ABSTRACT 

The  investigation  is  devoted  to  a theoretical  3tudy 
of  acoustic  wave  field::-  arising  frer.  a single  frequency  point 
source  in  slightly  viscous,  layered  liquids.  The  acoustic 
fields  are  derived  in  their  nodal  representation  for  three 
particular  sh. allow  water  eo ifiguraticn3i  one  ar.d  two-layex’ed 
media  bounded  by  perfectly  reflecting  planes,  and  the  two- 
layered  semi- infinite  medium.  The  resulting  node  foma  are 
compared  with  respect  to  orthogonality,  completeness,  finite- 
nos3,  medal  identity  and  discrefonoss,  physical  representations, 
and  "cut-off. ,!  Particular  attention  is  given  to  an  analysis 
and  extension  of  the  research  of  Br.  C.  L.  Pekeris  on  the  two- 
layered  semi-infinite  medium.  t derivation  of  the  general 
power  and  er.orgy  orthogonality  conditions,  based  on  studies 
of  electrceiagnctic  wave  modes,  is  also  presented. 
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the  propagation  of  explosive  sound  in  aha 1 loir  water.  In  this  paoer 
’V.  Pekeris  developed  a theory  of  acoustic  wave  propagation  in  a layered 


liquid  medium  and  successfully  applied  the  salient  features  of  the  results 

2 


to  some  experimental  data  obtained  by  Drs.  b'orzel  and  Swing.  In  hi3 
theoretical  work,  Pekeris  represented  the  fields,  arising  from  a single 
frequency  point  soui  e in  a horizontal  liquid  Layer  underlain  by  an 
unbounded  bottom  of  different  density  and  plane -wav 3 velocity,  by  a sum 
of  acoustic  "modes. " and  a .non-vanishing  "branch- line"  integral.  Although 
c.he  representation  was  adequate  for  his  practical  applications,  IV.  Pekeris 
left  unsolved  the  theoretical  problems  of  modal  orthogonality  and  justi- 
fication for  the  existence  of  the  branch  line  integral.  The  search  for 
the  answers  to  Hess  questions  led  to  a complete  re-examination  of  the 
entire  problem,  i'ho  results  are  embodied  in  this  paper:  a new  derivation 
of  the  field  integrals,  a n^w  transformation  and  representation  of  the 


Inal  solution,  u detailed  discussion  of  the  physical  interpretation  of 


bho  mathematical  results,  and  an  elaboration  of  some  general  properties 
of  acoustic  nudes,  in  tl.-o  muili-layereb  system. 

The  oovo I oo.mr: • t begins  in  fiyipter  2 with  the  derivation  of  the 
fundamental  acoustic  press'  . • vT  ^article  velocity  field  and  wave  equations 

for  a vertically  stratified,  slightly  viscous  medium.  Tho  snail  viscous 
effect  was  carried  through  only  as  a first  order  correction  tern.  A second 


^Pekeris  [ 12].  The  bracketed  numbers  refer  to  the  entries  in  the  reference  list, 
^t'orzel  [17]. 
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-joint  source  in  three  particular  media  are  derived  and  expressed  as  a aura 
of  modes.  The  first  two  media  are  one  and  two  layered  homogeneous  liquids, 
bounded  by  perfectly  reflecting  planes*  the  third  is  the  Pekeri3  config- 
uration. The  study  of  the  first  two  problems  provides  a convenient  means 
for  developing  relations  which  c m be  compared  with  the  Pekeris  results. 

The  examination  of  these  problems  also  suggests  the  form  of  tho  solutions 
to  be  expected  in  multi-layered  madia.  The  modoo  in  the  throe  systems 
are  compared  with  jesm^ct  ts  orth.ogone.lity,  completeness,  finiteness, 
nodal  identity  ar.d  disex-eo-.no.-.s,  physical  representations,  and  !l cut-off . " 

The  derivation  of  the  solution  is  given  at  the  beginning  of  each  of  the 
sections,  and  is  follows  1 by  a descriptive  discussion  of  tho  more  important 
features  of  the  solution. 

whereas  the  Investigation  in  Chapter  3 is  primarily  concerned  with 
the  modes  as  they  occur  in  particular  media,  Chaptor  4 is  devoted  to  the 
3tudy  of  soms  characteristics  of  acoustic  nodes  that  are  generally  valid , 
Particular  emphasis  is  plaeod  upon  the.  derivation  of  the  nouer  and  energy 
orthogonality  cor '?i tiers  and  the  physical  interpretations  of  the  propagation 
factois  from  energy  consideration?. 

Thoigh  tlv-  problems  are  restricted  t'  vertically  stratified  systems 
in  the  cylindrical  co-ordinate  system,  tho  methods  ran  be  readily  adapted 
me  general  v.v/o  problems  in  other  co-ordinate  systems.  as  long  as  the 
stratification  dees  not  lie  along  the  preferred  direction  of  wave  propaga- 
tion. The  general  source  distribution  as  u function  of  tiri9  and  space  can 
be  handled  by  well-known  courier  Integral  methods  combining  solutions  of 


i ' 

i\ 

it 


II! 

y-  it 

i £ f* 

¥~ 


• i mr«*M 


- X-S-J£«~-r-'-^£?  XlK+Xm  t 


the  single  frequency  point  source. 

The  application  of  the  results  given  in  this  paper  to  the  numerical 
calculation  of  acoustic  fields  in  shallow  water  problems  is  severely  limited 
by  the  idoalized  representation  of  the  medium  as  a system  of  homogeneous 
liquid  layers.  For  this  reason,  it  was  felt  that  numerical  examples  using 
average  oceanic  data  would  bo  superficial*  It  is  hoped,  however,  that  the 
discussion  will  fill  the  need  for  a relatively  simple  explanation  of  the 
physical  processes  occurring  in  shallow  water  sound  transmission  media  as 
represented  by  these  idealized  models. 

The  historical  background  for  this  research  is  contained,  for  the  most 

1 2 
part,  in  Pekeris's  article*  and  ir,  the  publication  of  Dr.  J.  lde,  et  al. 

Although  his  theoretical  work  has  been  entirely  replaced  by  Pekeris's  contri- 
bution, Ide's  descriptive  material  is  a good  introduction  to  the  problem. 

The  !IRC  Reports  are  readable  summaries  of  the  general  activity  In  submarine 

o • ^ 

acoustics  undertaken  by  the  United  States  Government  during  dor Id  War  II. 

Several  selected  references,  oach  having  an  extensive  bibliography, 

are  suggested  as  complementary  reading  material:  ?..  B.  aclcr*  on  the  subject 

of  inhomogeneous  electromagnetic  waveguides,  an  excellent  report  which  has 

greatly  aided  the  author:  H.  U.  M:.re'n,  Jr.  on  anomalous  eea-iic  wave  propa- 

' 7 

gabion,  M.I.T.  Radiation  Laboratory  Series,  Volume  15'  on  anomalous  radio 


^-Pekeris  [12],  p.  53. 

2Ide  [2j. 

s 

^Rational  Research:  Council  [3]. 

^National  Research  Council  [o]. 

5 Adler  [l] . 
r o 

f*cu.  Oil 

7M.I.T.  [5]. 
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CHAPTER  II 

THE  DERIVATION  OF  THE  FUNDAMENTAL  EQUATIONS 
The  orientation  of  the  cylindrical  co-ordinate  system  used  in  this 
investigation  is  shown  in  Figure  1.  Positive  % and  r are  denoted  the 
"depth"  and  "range"  respectively.  It  is  assumed  that  all  fields  are  inde- 
pendent of  the  azrauthal  co-ordinate  0. 
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The  acoustic  medium  is  assumed  to  be 
stratified  in  the  sense  that  its  char- 
acteristics are  functions  only  of  the 
F1c^  !•  depth. 

The  symbols  are  defined  as  follows: 
lo  Medium  parameters.,  time  and  frecuency  independent: 

Pq'z)  ^ Equilibrium  pressure, 

p (z)  ^ Equilibrium  density. 

c(z)  ^ Proportionality  factor  for  adiabatic  wave  motion. 

When  constant,  this  quantity  represents  the  velocity 
of  a plane  wave  in  the  unbounded  medium. 

F(z)  ^ Vector  body  force/unit  mass. 

ju(z)  ^ Viscosity  coefficient. 

2.  Field  quantities: 

*0*  P (z)  + InrjZ.t)  ^ Total  instantaneous  pressure,  where  'p'^ 

® time  dependent  excess  pressure. 

A A 

p1  = P (z)  + p(r,z)  * Time  independent  total  pressure,  where 

time  independent  excess  pressure. 

~p'  * p(.z)  + p(z)  s\'r,z,t)  ^ Total  instantaneous  density,  where  "s 

tine  dependent  condensation. 

p1  *■  p(z)  + p(z)  3(r,u)  ^ Tine  independent  total  density,  where 

s ^time  independent  condensation. 


v '..r,?.^)  = Instantaneous  vector  particle  velocity^ 
v (r,~)  ^ Tine  independent  vector  particle  velocity. 

The  source- fre^  "ovier-Stohes  equation  fo:  hydrodynamic  flow  is 

postulated  for  regions  in  which  the  nediur  parameters,  field  variables, 

2 

and  their  derivatives  are  continuous,  finite,  and  single- valued* 

.ay?  + /s'f-  vZ  =p'E'-Vp'+  4-/4V(v-2)  + *(v-V)£.  1.1 

C c>t  *•  -f 

The  following  assumptions  are  reasonable  for  acoustic  problems  of  interest 
here:^ 

1.  The  equilibrium  pressure  dependence  upon  depth  is  negligible. 


Therefore,  from  a consideration  of  tho  equilibrium  state,  pF  ~Vt£^Q. 

2.  The  pa.rticie  velocity  vector  is  irrotational.  This  is  equivalent  to 
neglecting  shear  waves,  an  assumption  which  is  reasonable  when  the 

viscous  terms  in  1.1  are  relatively  snail.  Thus, 

V x v sO  or  V (v*  y?)  — V'fV)2.)v 

3.  j5.5,v  are  snail  quantities.  Products  of  these  terms  or  their  derivatives 
arc  negligible  to  first  order.. 

Equation  I«1  on  the  basis  of  these  assumptions  reduces  to: 

Vp  = -fj£  1.2 

which  is  equivalent  to  teuton's  equation  of  motion  in  a lossy  (viscous) 
acoustic  medium.  The  equation  of  maso  continuity  for  a constant-mass  system 

i-For  simplicity,  the  particle  velocity  field  will  often  be  termed  the 

"velocity-field,"  and  mhe  adiabatic  fector  c(z)  the  "velocity  parameter. " 
'Page  [11],  Ch-vp.  S. 

Marsh:  [4],  p.  C. 


I 

» i 


it 


By  virtue  of  the  previous  asaumotions.  this  equation  ta'*es  the  reduced 


form: 

v p*  = ~pj§  • ^ 

A relationship  between  the  pressure  and  the  condensation  arises  from  the 
usual  assumption  that  the  wave  motion  occurs  adlabaticaliy,  which  implies 
that  the  excess  pressure  is  proportional  to  the  condensation: 

p = ec  s . 1.5 

The  condensation  is  eliminated  by  the  combination  of  equations  1.2,  1.4, 
and  1.5  whereupon  the  basic  acoustic  field  equations  result:- 


VP  = -P$f  +§/*?(V-:£)  , 1.6a 

V-P?  = • 1.6b 

czdC 

It  is  observed  that  the  viscosity  effect  is  found  in  the  force  equation 
but  does  not  enter  the  equation  of  continuity. 

The  generalized  "wave"  equations  are  found  by  performing  the  diver- 
gence operation  on  equation  L.6a  and  the  gradient  operation  on  equation  1.6b? 


v(v-pv)  =.  &.  <££  - _ ipfv±) 

c**dt*  3 c**  c)t\  C^J 


where  the  substitution: 


from  equation  1.6bywas  introduced.  The  solution  of  the  wave  equations  in 

Xsince  the  particle^  velocity  vector  is  irrotationa.l,  there  always  exists 
a velocity  potential  which  can  be  manipulated  as  a scalar  field. 

However,  for  the  problems  considered  here,  the  use  of  the  velocity  potential 
affords  little  added  convenience. 


8 


^T«ILV..dn»r-  ■"  « 


the  (general  case  Is  clearly  a difficult  analytical  problem  because  of  the 
complicated  interdependenco  of  the  field  variables.  But  for  the  medium 
consisting  of  stratified  homogeneous  liquids  the  wave  equations  are  simoly 
these i 


v’p~±<!£,  + £^>v3'4f-o , 

' 3 pc5- 

7¥-±£?-h±jL 

c*4tx  3 pcz  Zt 


1.7a 


1,7b 


where  the  medium  parameters  are  considered  constant  over  n given  Layer. 


.in* 


For  harmonic  tine  variation  C , the  field  and  wave  equations  read 


X7  _ 

*■  • -'C/ 


*** 


..3a 


pc*' 


v*P  + sd 


= o 


\ + ' Jj < = 0 


1,3b 

1.3c 

1.8d 


A 4 ait 

where  L = g — n , and  equation  1.7b  was  substituted  into  equation  1.6a. 

pcfc 

In  cases  of  practical  interest,  the  loss  factor  L is  snail  compared 


with  unity*  expansion  to  first  order  in  L provides  the  final  fora  of  the 
fundamental  equations: 


VyC?  =:  -ttUf>  (/~t  L)  J/ 


1.9a 


V’V_  =.  — LUi  p 

pc*r 


1.9b 


v%p  + a£(, 't-‘L)p 

V*JC  + “4(/-cL)* 


• O 

o . 


1.9c 


1.9d 


9. 


The  lo9 3 factor  thus  enters  the  field  equation  1.9a  as  a "complex" 
density,  and  the  wave  equations  1.9c  and  1.3d  as  a "complex"  velocity. 

The  complex  velocity  parameter  i.3  readily  interpreted  for  the  case  of  plane 
wave  propagation  in  the  infinite  homogeneous  medium.  The  plane-wave  expo- 
nential  representation  e has  a propagation  factor  .o(, 

whose  magnitudes  can  be  shoun  to  be: 

A * uJ  J 'j3-  f oC  — Co  ^ C / 

As  expected,  the  loss  factor  in  the  complex  velocity  parameter  gives  rise 
to  exponential  damping  and  a very  small  modification  in  the  free  space  wave 


number. 


The  boundary  conditions  of  continuity  for  the  wave  fields  are  specified 
from  purely  physical  considerations  as  follows: 

Consider  uwo  liquid  media  which  are  in  contact  at  the  equilibrium 
plane  * 0.  The  nediun  parameters  are  specified  as  shown  in  Figure  2. 

iO  r / 

< -i  j »*#.,  *-f 


Figure  2. 


2*0 


Across  any  surface  in  the  medium,  and,  in  particular,  across  the 
interface  ~ 0: 

A.  The  pressure  is  everywhere  continuous,  precluding  infinite 
acceleration  of  a finito  mass.  1.10 

3.  Except  at  a source  of  energy,  the  normal  component  of  particle 
velocity  ia  continuous  in  order  that  continuous  surface  contact  i3  insured. 

These  conditions  are  sufficient  for  the  determination  of  the  unique 
solution  of  equations  1.3.  No  specification  on  the  tangential  component  of 
particle  velocity  is  necessary  since  shear  waves  have  been  neglected. 


is  discontinuous  at  a discontinuity  of  the  complex  density.  This  statement 


is  interpreted  physically  as  the  fact  that  a moving  interfaco  can  support 
a pressure  field  which  is  continuous  but  may  change  dl3continuousiy.  On 
the  other  hand,  the  fundamental  mass  continuity  equation  1.3  is  not  defined 
over  the  interface  at  a density  discontinuity  3ir.ee  the  divergence  of  the 
discontinuous  density-velocity  product  does  not  exist,  even  in  the  presence 
of  viscous  effects.  This  essentially  means  that  the  mass  transferred 
across  the  interface  i3  not  conserved  since  the  moving  contact  surface 
involves  the  same  particle  displacement  of  different  masses  on  either  side 
of  the  boundary.  The  wave  equation*  therefore,  are  not  applicable  on  the 
interfacial  planes  between  liquids  of  non-uniform  equilibrium  density. 

These  peculiarities  account  for  the  rather  interesting  orthogonality  condi- 
tions to  be  discussed  in  Chapter  IV. 
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CHAPTER  III 


THE  SOLUTION  OF  PARTICULAR  PROBLEMS 

Having  established  the  fundamental  relations  among  the  field  variables 
and  the  medium  parameters,  we  now  concentrate  on  the  determination  of  the 
wave  fields  arising  from  a single  frequency  point  source  in  three  particular 
stratified  media  c Tlio  main  concern  is  the  derivation  and  analysis  of  the 
fields  in  the  acoustic  mode  representation.  The  investigation  lc  aimed 
toward  the  following  objectives : 

1.  To  extend  the  analysis  of  Or.  Pekeris  on  the  mode  solutions  in 
the  two- layered  unbounded  medium. 

9 T*  /3/»  » n 1 1 m •»  * V m>4ip  n O /w  nw  nl  /./«  f W a 

W (.  A xuc  U UUUU^XOa  OiiV.  V/Wmo^X»U  iJUUUJf  OA  lkU>4.ba  \J±  uwo 

three  major  forms  of  acoustic  modes  in  vertically  stratified  liquids. 

3.  To  demonstrate  the  power  of  the  Fourier-Bessel  integral  trans- 
formation method  of  solution  and  the  licited  application  of  the  well-known 
orthogonal  mode  expansion  technique. 

4.  To  present,  where  possiblo,  a reasonable  explanation  of  the 
phyaical  processes  involved. 

The  formal  problem  is  the  determination  of  the  pressure  and  velocity 
fields  due  to  a point  source  in  a layered  liquid  medium  when  the  fields  satisfy 
particular  conditions  at  tho  boundaries  and  at  the  source  location.  It  Is 
known  that  for  this  type  of  problem,  the  unique  solution  can  be  found  if 
either  the  pressure  of  the  normal  component  of  particle  velocity  is  specified 
over  the  bounding  planes.^-  This  paper  will  be  concerned  only  with  the 
solution  of  the  inhomogeneous  form  of  the  pressure  wave  equation  1.9c  and 
the  subsequent  determination  of  the  particle  velocity  field  by  the  use  of 
equation  l*9b. 

^Stratton  [16],  p.  485. 
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One  possibly  ambiguous  aspect  of  the  uniqueness  quality  of  the 
solution  is  this:  although  the  solution  is  uniquely  determined  by  eithe; 
the  pressure  or  velocity  designation  on  the  boundaries,  both  conditions 


I 

1.10  will  be  applied  at  non- totally- reflecting  surfaces.  The  difficulty 


is  resolved  by  noting  that  these  conditions  of  continuity  are  concerned 


with  the  total  resultant  field  at  the  boundary,  from  which  we  find  the 
fields  that  exist  independently  on  either  side  of  the  boundary.  Once  the 
field  is  found  everywhere,  the  uniqueness  theorem  states  that  there  is  no 
other  possible  solution. 

It  is  necessary  to  add  these  supplemental^'  conditions  on  the  wave 
fields  to  insure  uniqueness  when  the  diiim  extends  to  infinity*'*" 

A.  The  amplitude  of  the  total  field  approaches  zero  at  large 
distances  from  the  source  along  the  medium  co-ordinate  which  extends  to 
infinity.  2.  la 


B.  Active  acoustic  energy  that  is  transferred  must  be  diverging 
from  the  source  when  observed  at  large  distances  from  the  source  along  the 
medium  co-ordinate  which  extends  to  infinity.  2.1b 

The  Single  Homogeneous  Layer 

The  singlo- layered  medium  bounded  by  perfectly  reflecting  planes  is 
t’hown  schematically  in  Figure  3.  The  medium  represents  a shallow-water 
layer  bounded  by  a zero  pressure  surface  at  z » 0,  and  an  acoustically  "hard" 

^Sornorfeld  [14],  p.  188. 
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surface  at  z - H on  which  the-  vertical  component  of  particle  velocity  ia  0. 

The  layer  extends  to  infinity  in  the  radial  direction.  The  density,  velocity, 
and  viscosity  parameters  ore  constants  throughout  the  medium. 
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Figure  5. 
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A point  source  Q c “ which  has  a total  outflow  of  1 unit  volime/second 
is  located  at  (0,d).  The  source  as  it  appears  in  the  inhomogeneous  fora  of 
the  pressure  wave  equation  is  represented  by  the  product  of  Dirac  delta 
functions :~ 


where: 


Q ~ - KUJp(l’l  L)  '4)$(r) } 

/ (i  -a)  = oj  £?<*  ; f f(d  '£)*£  -X 
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2.2a 

2.2b 

2.2c 
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The  exact  statement  of  the  problem  is  the  determination  of  the 
solution  of  the  inhomogeneous  pressure  wave  equation: 

57/?  + “£(/- « *-)/*=■  “« * 2.3 

and  the  determination  of  the  velocity  field  fran  the  relation: 


Worse's  i’orra,  with  the  substitution  of  the  complex  density  [6],  p.  512. 
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subject  to  the  boundary  conditions} 


A. 

At 

3 75  0: 

p = 0 

2.4a 

3. 

At 

z - lit 

V §£  M 

0 

2.4b 

C. 

At 

large  r 

, the  total 

field  tends  to  zero. 

2.4c 

D. 

At 

large  r 

, the  field 

represents  divergent  radiation. 

2.4d 

if  radiation  exists.  (The  particle  velocity  in  component  fora  is 
Vz  » where  i0  and  £ are  unit  vectors  along  the  positive 

direction  of  the  co-ordinate  axe3.) 

Mo  shall  first  solve  this  problem  by  the  method  of  orthogonal  modes 
which  can  always  be  used  in  the  homogeneous  layer  configuration  when  the 

X 

homogeneous  conditions  are  satisfied  over  the  boundiag  planes. 
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2.5 


This  equation  represents  the  familiar  impedance  boundary  condition  which 
will  bo  briefly  discussed  later.  The  method  is  outlined  as  follows: 

1.  An  infinite  set  of  discreet  solutions  of  the  source-free  wqve  equation 
is  found,  all  members  of  the  set  individually  satisfying  the  boundary  condi- 
tions. 3acb  of  these  solutions  is  a ninodo.'' 

2.  These  modes  are  shown  to  comprise  an  orthogonal  set  of  functions  whose 
properties  provide  a representation  of  the  source  term  by  a linear  combina- 
tion of  the  modes. 

3.  The  pressure  field  is  then  found  to  be  another  linear  combination  of 
these  modes,  whose  amplitudes  are  evaluated  in  terms  of  the  source  expansion 


^Somnerfeld  [14],  p. 
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coefficients 


In  cylindrical 

<?£  X 
JrK  T r 


co-ordinates,  the  homogeneous  form  of  equation  2.3  ie: 

+ 2.s 

t>Z  C 4 


If  it  is  assumed  that  the  presatire  field  can  be  expressed  in  the  separated 
form: 


p(r,  = Rlr)F(z)  , 


then  equation  2.6  is  directly  reduced  to  two  ordinary  differential  equations 


Sr  + j.iS  . FR  - o 

<£F  t (sU-lJ-  Fjf  ~o , 
T (c*  y ; 


2.8  a 


2„3b 


,3 


whore  % is  the  separation  constant.  For  later  convenience,  let: 

k + id,  k and  a real,  If™  Propagation  factor. 

'p  - Distribution  factor. 

The  s co-ordinato  boundary  conditions  ar«  ncu  applied  tc  2,3b.  The 
result  is  an  infinite  number  of  diocre-l  solutions  of  2.6,  defined  by  the 


relations: 


F = RJr)Fh(i)  , n * />', 
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Now  let  us  assume  that  the  z dependent  delta  function  can  be  expanded 


into  or.  infinite  series  of 


The  constants  A can  be  evaluated  by  application  of  the  orthogonal  properties 


of  the  F- functions 
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2.11 


where  the  fiolds  in  the  range  Z/*H  are  everywhere  0.  An  integration 


of  equation  2.10  over  the  range  of  z determines  the  constants: 

A 


' 7?J  *' — 


This  series  representation  of  the  delta  function,  however,  does  not  converge 


for  any  value  of  z.  The  justification  for  its  use  lies  in  the  fact  that  the 


final  representation  of  the  field  as  a series  does  converge  for  all  r and  z 


except  at  the  source  singularity  r = 0.J 


inn  3oxuvxon  oi  ni.e  xnnamogexieous  wave  equation  £.6  is  now  assumed 


to  bo  a linear  combination  of  modes i 


ptrt-S-tyrifJi), 


2.  .13 


in  which  each  of  the  functions  R contains  an  unknown  constant  factor  which 

n 


will  be  evaluated  in  terms  of  the  source  expansion  coefficients.  Introducing 


thi3  expansion  and  the  preceding  relations  into  the  wave  equation  2.3,  we 


can  produce  the  following  equations: 


-21  Fn(. 2)( ^ 1 


2.14 


A similar  situation  exists  in  on--  v.ibrati rr 
Morse  [c]  n.  93. 
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string  problem  discussed  by 


The  application  of  the  orthogonality  conditions  2.11  provides  t.ho  equations 
Tor  the  R- functions: 


0*  t «X 


let:  0 LJ)  Fn  (*t)Fn  (r)  . 

Then,  jyi^  + - 


1.15 


The  Hn  functions  are  thte  Green's  functions  fcr  the  cylindrical  system  which 
represent  diverging  radiation  and  vanishing  amplitude  at  large  ranges*  and 
satisfy  the  singularity  requirement  at  r 3 9.^  Corre spun ding  to  the  time 
""actor  the  aooronriate  solutions  of  2.16  are  t.hn  zero  order  Hanlcel 

Functions  of  the  second  kind  with  a normalizing  factor  of  Ti.  That  this 
Hankel  function  represents  diverging  radiation  is  3een  from  its  asymptotic 
form  at.  large  ranges: 


>%r 


2.16 


The  pressure  and  velocity  fields  are  now  completely  determined: 

P ~ SmpJ  s,'npnir  KHr)  2.17a 

X * ;"pn  2 Ho'fth  >*))  . 2.17b 


A discussion  of  the  solution  will  be  undertaken  later. 

Ve  shall  now  3olvo  the  preceding  problem  by  the  Fourier-Bessel 
integral  method  which  alone  will  be  applicable  in  the  later  problems.  This 
method  is  considerably  more  powerful  than  the  method  of  orthogonal  modes 
since  medal  orthogonality,  among  other  qualities  of  the  modal  set,  is  not  a 

florae  [7]. 
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prerequisite  for  a solution  but  rather  a consequence.  The  technique  in 
outline  form  is  as  follows: 

1.  Too  point  source  hem  is  '"pressed  us  a Fourier-Besael  transform 
which  physically  represents  a linear  combination  of  a continuum  of  plane 


sources. 


2.  The  pressure  field  due  to  the  presence  of  one  plane  source  is 


found. 


3.  The  field  of  the  infinite  continuum  of  plane  sources  acting 
simultaneously  are  lir.oarly  combined  to  give  the  field  of  the  point  source. 

Any  reasonable  function  of  the  range  co-ordinate  can  be  expressed  in 
an  integral  form: 


/( <)  (kr)k  *tk 


when  the  integral  existsT  J^kr)  is  the  zero-order  Bessel  function  of  the 
first  kind,  k i3  purely  real,  and  g(k),  refined  as  the  Fourier-Bessel  trans- 
form, is  found  from  the  inverse  transformation: 


Let  us  express  the  radial  delta  function  as  a Fourier-3essel  integral.  Its 
transform  is: 

oO 

f(k)  — ft 'r)dl (ktjtJf  - 1 2.20 

and  the-  final  expression  for  the  Fouricr-Bessel  integral  is: 


J(r)  — f-Z  (kr)  k Jk 


It  is  easy  to  shew  from  an  integration  by  parts  that  the  above  integral  does. 


^Socrnerfald  [14],  p.  £40. 
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in  fact,  diverge  for  all  r.  This  representation,  however,  will  be  used  by 
analogy  with  the  previous  case  of  the  non-converging  series  representation 
of  the  z delta  function,  where  the  justification  is,  again,  that  the  final 
expression  is  everywhere  convergent  except  at  the  source  axis,  r ® 0. 

The  source  tern  of  equation  2.5  in  its  integral  fom  la: 

£X? 

Q s=  (~~LU>p(l  “<  dk.  , 2.22 

o 

V.’e  can  give  a direct  physical  interpretation  to  the  above  expression.  The 
integrand  represents  a plane  acoustic  sourco  at  z c d,  with  a radial  amplitude 
distribution  of  magnitude  TJ  (kr).  Figure  4 is  a schematic  representation 
of  the  plane  sourco  distribution. 
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bet  us  now  assume  that  u single  plane  source  exists  in  the  medium  at 
the  phone  z a d with  a radial  amplitude  distribution  corresponding  to  the 
number  k.  '-.’e  shall  first  find  the  fields  due  to  this  source  by  solving  the 
wave  equation: 


+ * -lcu^( PlL) %( i -cj) A J^(kt')3 


2.23 


subject  to  the  sane  boundary  conditions  as  before,  equations  2.4?-  The 
oressuro  field  Is  assumed  to  be  composed  of  two  temst^£  ^ is  the 


^This  method  v;as  suggested  by  the  derivation  given  by  Adler  for  the  case  of 
electromagnetic  radiation  from  a dipole  in  a cylindrical  dielectric  rod 
which  exists  if  free  space  LlJ*  App.  4. 
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''particular11  solution  of  the  wave  equation!  it  completely  represents  the 
presence  cf  the  source,  and  only  tho  source,  and  is  chosen  a3  though  the 
medium  were  unbounded.  i3  the  "complementary"  solution  of  the  source  free 
wave  equation  and  13  so  combined  with  ^ that  the  total  field  satisfies 
the  boundary  conditions  at  the  reflecting  plano3.  It  is  explicitly  atat-sd 
that  neither  of  these  terras  satisfies  the  boundary  conditions  independently. 

Tho  field  transition  conditions  across  the  source  plane  rau3t  first 
be  established.  Lot  the  fields  above  and  below  the  plane  z •*  d be  denoted  by 
the  subscripts  1 .and  2 respectively.  Since  the  pressure  is  everywhere 
continuous, 

Af  2 - J : A.  z-  . 2.24 

' ~ Tp,  fpJL 

However.  tho  particle  velocity  vector  is  not  continuous  because  the  fluid  is 
moving  in  opposite  directions  above  and  below  the  source.  The  magnitude  of 
the  velocity  discontinuity  is  found  from  the  fact  that  the  vertical  component 
of  particle  velocity  at  the  source  surfacp  Is  equal  to  the  total  volume  out- 
flow/second  in  the  z direction.  Since  the  delta  function  represents  a volume 
outflow/unit  volune/second,  the  total  outflow  in  tho  z direction  is  simply: 


J 

o 


f^f{i'J)kJ0(/<r)Jt  - £<Z(Ar)  . 


p Or 
C • *-o 


The  normal  component  of  particle  velocity  at  the  30urco  surface,  therefore, 
has  the  magnitude: 

IvJ  =.  /cUl(kr)  . 

Consequently,  the  velocity  discontinuity  at  tho  source  plan9  i3: 

At  z =J:  - bul  - -ZkZ(kr)  } 2...K 


^The  minus  sign  is  chosen  for  correspondence  with  the  sign  of  the  downward 
traveling  pros jure  wave  in  equation  2.27b. 
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The  particular  solution  ^ can  now  ho  found  by  aasminj  anprooriate 
solution-  of  tV  wave  equation  2,24  in  the  re  -ions  c.bove  and  bo  lev  the  source* 
qnu  "latching  the  fields  at  t-V  piano  2 3 d.  ..‘inoe  the  boundary  conditions 
will  be  applied  to  the  total  field,  wo  run  choc 30  a solution  which  represents 
upward  traveling  waves  above  the  source  arid  downward  traveling  waves  bo low 
the  source.  The  acccotable  choice  is: 


^4  = 


provided  eh-  1. 


— A e L-^?r  T>(k r)  / Ori it  ~ J 

= 8e  Cp(z  ~J)TAk ?),  J*  . 

3?  - toZ{/  ~lL)  -k'2'- 

r ~Ta.  ( 


2.27a 

2.27b 

2.28 


The  constants  A and  D arc  evaluate-;  by  "lutchi.nj  tun  above  equations  ever  the 
olan9  7,  ~ d by  the  application  of  equations  2.24  and  2.26.  The  particular 
solution  is  found  to  be: 


= - 1 uJf  (/- 1 L)^f^  U %(&r)  , 0£  £■  — d/ 
» - * fj'i L) £ <X (A y) > 26//, 


> On 
!s-  p 


The  conp loitnutary  solution  is  sinply  c honor  as 


/ = Me  ? gJ0'(A  rj  rNe^ZM, 


vhicn  represents,  respectively,  upward  and  downward  tr.velirr;  waves  in  the 
absence  of  a source. 


The  constants  are  Co  end  iivoi.  ily  fro.  i the  application  of  the  s-co-ordinat 3 


boundary  conditions  2.4a  and  t to  tho  sur.  of  'J-'p  and  * 


M - Dke  ^ 5(0/3 J ■ /V  - -Ok  cos  pf/i-tf)  t 


2.51a 


where 


&C0S/3//  ' 

jO  it  - t (0/0  ((  - t L ) 


2.31b 
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The  conditions  2.4c  and  d will  be  applied  only  to  the  final  solution 
since  tuey  are  compatible  only  with  source  distributions  located  entirely 
within  a Tiniie  radius  of  the  origin.^- 

The  fields  of  the  criminal  point  source  are  found  by  integrating 
equation  2.13  over  the  complete  range  of  k»  that  ia  to  say,  we  aun  the  pressure 
fields  due  to  the  contributions  of  an  infinite  continuum  of  plane  sources 
operating  simultaneously  over  the  plane  2 * d: 

M?  0t> 

+ «£(/ -'Qjjjtst/c  » -in/p(/-ti)£(2-J)f  fr/,  2.32 

o ^ '■» 

The  field  of  the  point  source  in  tho  single-layered  medium  is,  then: 

. «p3 

p - i%(kr)Jk  + fe 
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J0  piospH 
t ; o i=  £±4 


The  first  integral  must  obviously  represent  the  source  completely  and 

3hould  be  expected  to  reduce  to  th^  field  of  a point  source  in  the  unbounded 

medium.  That  this  act-ua'7,.  vseurs  is  shown  by  an  extension  of  the  Sconierfeld 

o 

transformation  to  tho  lossy  case:1"  /, 


aon  00  one  xossy  cu.3e:  /, 

**  , , /uj*r,  •/  if*.  j^lr* 

V '"-few  - . J ■ 'i 


2,54 


Cr*r(?-4)*J& 


For  ease  in  rianipuiaiion.  the  expressions  for  tho  totu-  field  con  be  contracted 
into  the  following  fores: 


^'3t  ration  ^16],  p.  485. 
^lonnerfeld  j 14] , p.  240 
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p -J' lDTAkr)k<lk.  COSB/fccjp  J 4- sjnpftstnpjj, 
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p ~f  2£>T*(fcr')k<fk  iihpJfcozpHctSpZ  rl'OpH £tnB2) f 2.3t 

j 3 (,ospfT  ' cilFBtfi 


or,  in  tho  abbreviated  fashions 
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p'jJ.6<g)Tjkr)KJk 


2.36 


The  remainder  of  the  derivation  is  devoted  to  the  exact  interpretation 
and  evaluation  of  the  field  integrals  2.35.  Considered  as  functions  of  the 
real  variable  k,  the  integrands  are  not  analytic  everywhere  due  to  the 
presence  of  branch  points  and  simple  pole3.  V/e  shall  therefore  extend  the 
integrand,  by  means  of  the  orincinie  of  analytic  continuation,  into  the  complex 
plane  where  the  integral  can  be  properly  defined  by  an  appropriate  choice  of 
a path  in  the  complex  k plane.  To  this  end,  tho  conn lex  variable  ^ = k + ia, 
k and  a real,  is  substituted  for  the  real  variable  k. 

The  major  objective  is  to  transform  the  existing  integrals  into  line 
integrals  along  closed  pati.-s  about  tbs  poles  of  the  integrand  in  such  a 
manner  that  only  integrals  along  tho  polo  contours  contribute  non- vanish... sg 
terms.  Particular  attention  is  given  to  circular  contours  with  large  radii 
on  which  the  integrands  provide  vanishing  contributions  as  the  radii  approach 
infinity.  In  their  present  formulation,  the  integrals  will  nowhere  vanish 
for  largo  $ because  the  Bessel  function  grows  in  magnitude  for  large  values 
of  a complex  argument.  The  transformation 


iz(tr)  =/£%-;  t H~pr) 


Mb 


is  introduced  to  insure  that  tho  integrands  -.’ill  vanish  for  large  From 


their  asymptotic  forms, 
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it  it;  apou’-er.t  th  it  the  Ilar.l  el  function  of  t'-e  'Irst  •'inf  will  vanish  on  an 
infir.it*’  contour  in  the  inner  half  ^-plune,  and  the  liunkel  function  of  the 
second  kind  '..’ill  vanish  on  an  infinite  contour  in  the  lover  half  J-pInne. 

'•a  shall  next  clear  the  natter  of  the  doubLc- value -'ness  of  the  inte- 
grands of  2.25  duo  to  the  double  choice  of  the  sign  of  as  seer,  from  its 
definition  in  equation  2.2S.  bach  choice  of  sign  defines,  essentially,  ono 
branch  of  a two- sheeted  aeinunn  surface."*'  i'k*?  sign  of  the  root,  and  hence 
the  choice  of  sheet  on  which  the  integration  oath  is  to  iiov  is  chosen  so 


that  the  integrals  2.55  will  converge  for  l.-r-ge  k by  insuring  that  the  inte- 

2 


grand:;  rapidly  approach  aero-  as  k approaches  infinity.  ' occept  near  tine 
orsnch  point  P"r0,  it  is  observe*3  that  ® is  e.iso.itialiy  imaginary  for  large  k. 
The  positive  sign  of  the  rcct  is  thus  indicated  with  i dofintsd  as  follows: 


For  = fjsfO-iL]  - arj p ko 

A?'**.  ■ 3 c'fAz- zlU-‘LJ 3* -IF 

C ■ I c*-  /j  J, 


2.59 


holing  -the  restrictions  on  the  ranges  of  z.  me  can  easily  show  that  the 
jritegrunds  of  2.25  will  assuredly  vanish  for  largo  k. 

It  is  no:?  necessary  to  incart  a “branch"  cut  into  the  complex  plane 
in  order  that  the  anulyt icily  of  the  integrand  is  preserved.  The  cut  starts 
t the  branch  point  and  extends  to  infinity i as  long  as  the  phase  specifica- 
tions on  j-  in  equations  2.59  are  adhered  to,  the  cut  can  be  oriented  in  any 
direction.  A tentative  choice  in  a cut  parallel  to  the  imaginary  axis  as 
shown  in  Figure  5,  where  the  snail  imaginary  part  of  p-  has  been  neglected  in 
:he  drawing.1'’  Approximate  phases  of  J:  are  ;ivon  at  four  points  near  the  cut. 


morse  [7_t 

"Sorr.orfeld  [ 14] , p.  251. 

This  particular  cut  is  suggested  by  Pe’oris'n  method  [12],  Tart  2. 
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Regardless  of  the  branch  cut  orientation,  the  argument  of  {?  changes  190° 
across  the  cut-  This  is  easiest  seen  by  expressing  f in  the  form: 

3 - (*£  - *Z)'*‘  =T  ( C - _C , 

V ~ i<p  ~ ” O-'L)**- 

Defining  the  complex  r ector  ^ and  ^ real,  a -J4-/  f ve  observe  that: 

/✓  ; $ / !/  £ 
r>  l ^ o / . v 

B - h <?  -T  oj  , 

/ ' c 

It  ’.s  clear  that  ^ changes  by  2u  radians  for  choices  of  £j  on  either  side 
and  adjacent  to  the  branch  cut« 


The  next  operation  is  the  datamation  of  the  poles  of  the  integrands 
of  £.35  which  lie  at  the  values  of  Z for  which* 

CcS&/i  *0  rt  - J,Z  -'•£>  . 

The  poles  are  those- 
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~ k tl<-kr{<r>  " J£L0"1  L)  ~/3*  . 

CK 
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All  of  the  poles  in  the  right  half  of  the  plane  lie  in  the  fourth  quadrant 
in  th.9  region  k < u>/ c.  Fran  the  imaginary  part  of  the  above  equation  it  is 
seen  that  the  polos  lie  on  a hyperbola  which  approaches  Infinity  along  the 
negative  imaginary  axis.  In  the  lossless  case,  there  exist  a finite  number 
of  poles  cn  the  real  axis,  and  an  infinite  nuater  on  the  negative  imaginary 
axi  s. 
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The  field  integrals  will  now  be  defined  ir.  the  conolex  tf-olane  of 
Figure  5.  Introducing  the  Hankel  function  transformation,  2.37 ; the  integrals 
in  2.3o  take  the  form: 

(X> 

/ = [£(/>) H,  fa* Jk 

^ W 

Cauchy's  Theorem  provides  the  following  representation  of  the  above  equations: 

f& (a) Ho^ib) M/t  - f &(s)&J)(6r)¥4X  > 

6 JAB  i-BC 

•o 

I &(p) *-C(kr)k«k.  ^-  (& r)  VcM 

* JfiB+BF+FG+G-H  + Hrt  Poles' 

in  which  the  direction  along  the  path  is  given  by  the  order  of  the  subscript 
letters,  and  the  last  integral  represents  the  sum  of  the  integrals  on  contours 
encircling  ill  the  poles  in  the  fourth  quadrant,  taken  in  the  clockwise 
direction.  The  paths  CB,  FE,  and  Jli  are  chosen  to  he  segments  of  a circle 
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vhoa^  center  is  at  the  tulg in.  is  the  radius  of  the  circle  tends  to  infinity, 
the  integrands  along  the  paths  BC  and  HJ  ocor.e  vanishingly  small  for  ail  z. 
and  nil  r f 0.  l.’.onc  FF  however*  the  integrand  will  vanish  only  when  the 
range  is  at  least  of  order  of  the  depth  and  of  the  wavelength  since  the 

Imaginary  part  of  “ is  positive?"  This  restriction  is  rather  awkward*  a 

3 

ropresentati in  which  is  everywhere  valid  is  more  desirable.  This  aim  can 
be  achieved  by  choosing  the  branch-cut  orientation  so  that  the  imaginary 

r 

part  of  p is  negative  ever  the  entire  right  half  tf-planat  The  remainder 
of  this  paragraph  Is  restricted  to  the  lossless  case.  It  is  easily  shown 
fro  ; equations  1.19  chat  the  branch  cut  should  lie  along  the  path  fm(p)  = 0* 
more  precisely,  the  branch  cut  i?  along-  the  real  axis  from  0 to  _/c . and  along 
the  negative  iragiju.r.y  arils  from  0 to  oo.  r.s  shown  ir.  Figure  6. 
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The  phases  of  p arc  also  indicated  on  this  figure.  The  pole3  are  considered 
to  lie  -n  infinitesnal  distance  away  from  the  branch  cuts  by  the  assumption 
that  an  infinitesimal  loos  io  present,  however  slight.  Cauchy's  Theorem 
then  provides  this  representation  of  the  integrals  of  equations  2.35,  refer- 
ring the  contour  designations  to  Figure  6: 

oO 


The  paths  C8  and  FE  are  again  choson  to  be  segments  of  a circle-  whose  center 
is  at  the  origin,  and  the  integrands  along  these  paths  tend  to  zero  as  the 
radius  approached  infinity  for  all  z and  all  r f 0.  The  3um  of  the  integrals 
along  AC  and  GF  are  shewn  to  van!  ’h  by  Pekeris.^  The  integrals  along  the 
branch  line* 


since  p simply  changes  sign  on  either  side  of  the  cut  and  the  integrand  is 
even  with  respect  tc  p.  The  only  remaining  contribution  to  the  field  ccroes 
from  the  line  integrals  around  the  pole3,  which,  of  course,  is  the  sum  of 
the  residues  of  the  integral!  thus 


® zJT^  Res 


The  residue  at  the  nth  pole  % i3  given  by  the  application  of  a well-known 
theorem: 


2.41 


^Teksris  [12],  Part  2. 


The  residues  can  be  evaluated  easily,  and  the  result  1 the  same  for  fields 
both  above  and  belou  the  source: 

P - Mi  ^J)  H Si»pJsthj3„  a A vfoAr).  2.42 

H /)*/ 

This  expression  is  precisely  that  derived  by  the  method  of  orthogonal  modes, 
equation  2.17a.  The  particle  velocity  field  is,  of  course,  also  the  same. 
This  correspondence  atte3t3  to  the  validity  of  the  Fouri9r-3essei  method. 

Tie  effect  of  a small  loss  doe3  not  change  the  derivation  signif- 
icantly. A.  branch  cut  can  always  be  found  such  that  the  Imaginary  part  of  f 
is  negative  over  the  right  half  ^-planei  it  takes  the  form  of  a segment  of 
a hyperbola  starting  at  the  branchpoint,  and  extending  to  infinity  along 
the  negative  ir aginary  axi3.  The  line  integrals  along  all  the  cor '.ours 
except  th?  pole  contours  will  vanish  for  exactly  the  same  reasons  as  outlined 
above.  Tie  residues  have  the  same  form  as  equation  2.42,  except  now  the 
poles  are  all  complex. 

Toe  solution  will  now  be  discussed  ia  detail. 

1.  A comoarison  of  the  two  methods  of  solution  indicates  that  those  charac- 
teristics of  the  modal  field  representation  which  are  necessary  tc  success- 
fully apply  the  method  of  orthogonal  modes  are  actually  consequences  of  the 
solution  -/hen  derived  by  the  Fourier-Bessel  integral  method.  These  charac- 
teristics are  completeness,  orthogonality,  and  satisfaction  of  boundary 
conditions^ 

A.  Completeness  of  a modal  set  essentially  provides  the  guarantee 

^Suggested  practical  references  for  these  concepts  are  the  investigations  of 
Somrierfolct  [14]  and  Korse  [3]. 
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that  an  arbitrary,  continuous  acoustic  field  in  the  medium  can  be  exactly 
represented  by  a linear  combination  of  the  nodes.  The  set  of  nodes  was 
assumed  to  bo  complete  in  the  orthogonal  node  method,  whereas  completeness 
of  the  modal  set  was  insured  from  the  integral  transformation  method  in  view 
of  the  disappearance  of  all  cc  itour  integrals  but  the  residue  components. 
lie  shall  later  find  that  the  sot  of  discrete  aodes  is  not  complete  when  the 
medium  is  unbounded  in  the  z direction. 

2.  ’.lodai  orthogonal ity  is,  of  course,  the  major  prerequisite  for 
the  application  of  the  ortl  ogonal  mode  technique.  This  property,  however, 
is  clearly  a by-product  of  the  solution  by  the  integral  transformation. 

C.  A further  requirement  for  the  use  of  the  orthogonal  mode  method 
is  the  necessity  that  each  of  the  modes  satisfies  the  boundary  conditions 
independently.  In  the  Fourier-Besael  method,  the  initial  fora  of  the  solu- 
tion contained  two  toms  which,  did  not  independently  satisfy  the  boundary 
conditions  but  eventually  were  transformed  into  the  mo^al  set.  It  is  thus 
reassuring  but  not  requisite  that  the  integral  transformation  should  result 
in  an  infinite  3ot  of  discrete  solutions  of  the  wave  equation.  This  point 
is  quite  signif leant  for  solutions  in  medium  configurations  which  are  un- 
bounded ir:  the  7,  direction. 

D.  As  long  is  the  nediutu  is  homogeneous,  and  subject  only  to  tha 
homogeneous  boundary  conditions,  we  can  say  that  orthogonality  and  completeness 
Imply  one  another.  However,  in  the  multi-luyer,  or  inhomogeneous  case,  we 
must  examine  the  final  result:;  of  the  integral  transformation  method  in  the 
individual  case  before  a judgment  is  made.  La  the  next  two  problems  we  shall 
find,  respectively,  examples  of'  completeness  without  orthog  nality,  and 
orthogonality  without  completeness. 

2.  It  is  noteworthy  that  tt-  final  solution  determined  by  the  Fourier-Bessel 


•vMi  jjggyw? 


integral  transfomvitlon  is  essenti  ally  independent  of  the  choice  of  branch 


cut  since  the  pole  contours  alone  contributed  to  the  solution,  './hen  the 


medium  is  unbounded  in  the  7 directions,  «e  shall  pay  this  severe  penalty 


x or  extend  ing  the  integral  into  th3  complex  plane:  the  field  representation 


will  depend  upon  the  choice  of  branch  cut,  resulting  ir.  an  infinite  number 


of  available  representations. 


5.  Several  physical  interpretations  of  the  individual  modes  are  possible 


when  the  Kankel  function  is  approximated  by  it3  asymptotic  form: 


$/ni3£  e V 

' f , . \ ‘a  ..’A 


It,  is  clear  that  the  mode  represents  an  inhomogeneous  [z  dependent)  cylin- 


drical wave,  which  is  diverging  Iron  the  source.  This  wave  propagates  with 


the  chase  velocity 


y - — and  damps  exponentially  at  the  rate  of  nepers. 


unit  range,  in  addition  to  the  geometrical  spreading  indicated  by  the  factor 


r ' . In  the  lossless  medium,  thoro  ire  a finite  number  of  modes  for  which 


°^n  -3.  A term,  of  this  nature  represents  an  undamped  cylindrical  wave 


which,  as  3hown  from,  considerations  in  Chapter  IV,  radiates  energy  without 


loss.  A term  for  which  fin  w 0 represents  an  oscillatory,  damped  field  which. 


is  easily  shown  to  carry  no  real  energy. 


At  large  ranges  where  the  factor  f '*•  can  be  neglected  when  compared 


to  the  exponential  factor,  the  undamped  01  freely  propagating  mods  can  be 


interpreted  as  tbs  resultant  of  two  constructively  interfering  traveling 


"plane"  waves.  For  a given  mode,  the  plane  waves  are  traveling  at  discrete 


angles  @ with  resooct  to  the  vertical  given  by  the  well-known  relation; 


57  n Sn  ar  k„  c. 
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We  thus  f ind  that  the  plane  sources  with  propagation  factors 
operate  into  a "resonant"  condition  in  uhich  energy  can  be  transmitted  without 

1 *>+.*+*%  I’t,*,  — —— ._x  _ /»f* x _ r*  _ i *»  _ a.  i — _ i . . i . ji  . _ a.  __  _ , _ a.  i„._  p x . _ 
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ference  with  no  resultant-  radiated  energy. 

4.  It  uas  found  that  second  quadrant  poles,  representing  converging  waves 
did  not  appear  in  the  solution  because  the  contours  did  not  leave  the  right 
half  $-pione.  There  were  no  first  or  third  quadrant  poles  of  the  field 
integrands.  These  poles,  which  correspond  to  waves  whose  amplitudes  increase 
exponentially  in  the  radial  direction  are  excluded  as  a direct  consequence 

of  the  fundamental  equations  and  not  a3  a result  of  a separately  imposed 
boundary  condition.  In  Chapter  IV,  this  statement,  is  proved  to  be  applicable 
in  the  multi- layered  system. 

5.  It  is  recalled  that,  even  in  the  gonerally  lossy  case,  the  distribution 
parameter  is  real  and  a function  only  of  the  layer  depth  H.  The  modal 
dependence  upon  the  vortical  co-ordinate  is  consequently  a purely  sinusoidal 
function  which  is  fixed  when  the  layer  depth  is  specified.  It  is  therefore 
possible  to  use  the  and  the  resulting  vertical  field  distribution,  as 

a label  identifying  the  nth  mode.  It  is  particularly  important  to  note  that 
the  codes  exe  clearly  identified  when  considered  a function  of  source  frequency. 

6.  Let  us  again  study  the  lossless  case.  From  equation  2„40,  it  is  observed 

that  ail  of  the  poles  on  the  real  axis  fall  in  the  range  O*-  ~r~  » for 

all  frequencies.  It  is  also  clesr  that  kQ  decreases  from  at  the  high 
frequencies  to  0 at  a certain  finite  frequency,  which  is  defined  as  the 
"cut-off"  frequency  of  the  nth  mode.  Belov?  the  cut-off  frequency,  the  pole 
on  the  real  axis  disappears,  and  a now  pole,  corresponding  to  the  same  pn, 
appears  on  the  imaginary  axis;  the  node  thus  changes  from  free  propagation 


node  series  are  related  by  a Poisson  summation.  "Image- type"  representations 


of  fields  in  the  nulti-layered  system  are  possible  but  will  not  be  discussed 


in  this  paper,  since  Dr.  Pekerls  has  examined  them  in  detail. 


9.  The  solution  for  the  modes  in  the  system  with  the  general  homogeneous 


or  "irpedance"  boundary  conditic  n has  the  same  form  as  the  solution  given 


in  equ&tior  2.42.  This  condition  represents  a fixed  relation  between  the 


pressure  and  particle  velocity  fields  over  the  bounding  plane,  and  conse- 


quently, can  exist  only  when  the  wave  phase  fronts  are  approximately  plane 


or  where  rray"  theory  applies.  Dr.  Pekeris  has  derived  an  expression  for 


the  Unit  of  applicability  of  the  Impedance  boundary  aporoximation  In  the 


layered  Liquid  medium. *" 


The  Two- Layered  Medium  Pounded  by 


Reflecting 


The  configuration  for  the  bounded  two- layer  medium  is  shown  in 


Figure  7.  The  upper  layer  represents  a shallow  water  layer  bounded  by  a 


zero  pressure  surface  at  z * 0.  The  .*  ov;er  layer  represents  an  idealization 


of  the  ocean  bottom  in  which  it  is  assumed  to  be  a liquid  with  appropriate 


acoustic  constants.  Tho  liquids  are  in  continuous  contact  at  the  interface 


z * H.  The  lower  layer  is  bounded  by  on  acoustically  "hard"  surface  at 


z 31  H + h.  A point  source  of  energy  is  located  at  the  point  (0,d). 


Figure  7. 


^"Pekeris  [13]. 
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Tho  constants  of  the  lower  layer  are  Gsstmed  to  be  greater  than  the  corre- 
sponding constants  of  the  unppr  layer.  Tho  subscripts  1 and  2 refer  to 
quantities  in  the  upper  and  lower  layers  respectively. 

The  exact  statement  of  the  problem  is  the  search  for  the  solution  of 
the  pressure  wave  equations  in  each  of  thr  layers* 

(7/7  L,l  ffr)  O it  2.45a 

c< 

vh 

U*  2 ,45b 

and  the  determination  of  the  particle  velocity  fields: 


U 

¥ 

1 

n 

2.46a 

ssA  = vga. , 

2.46b 

subject  to  tho  boundary  conditions: 

A. 

.\fc  z " 0:  p^=  0. 

2.47a 

B. 

At  z " K:  py  and  v?1*  v^g. 

2.47b 

C. 

At  z 11  H + h:  v n3  0. 

2.47c 

D. 

At  large  r,  the  fxelds  tend  to  zero. 

2.47d 

E. 

At  largo  r,  the  fields  represent  diverging  radiation,  if 
radiation  exists. 

2.47e 

V/o  shall  now  show  that  the  method  of  orthogonal  mod03  is  not  applicable 
in  this  problem.  Following  the  established  procedure,  we  assume  that  the 
pressure  fields  have  the  separated  forms* 

fi  = /?,  tr)  fw  J ft  = /?  (r)£  U). 

The  homogeneous  forms  of  wave  equations  £.45  are  thereupon  reduced  to  ordinary 
differential  equations  of  tho  form* 
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provided  that: 


Ps  = £y(t  -‘Lsi  -*s  • 


2.43 


Tho  F-funetions  are  f irst  chosen  to  be  solutions  which  fit  the  specifications 
at  j - 0 and  H + h* 

fx  * £^/3^fz-[Htk3): 

the  unknown  constants  vi.1.3  be  evaluated  fron  the  boundary  conditions  at  n -di. 
It  is  sufficient  to  assume  that  a constant  factor  31103  in  the  R-functioos, 
and  take  the  constant  A to  be  unity.  The  continuity  conditions  at  H are 
satisfied  ulien: 


P - sia&.tLj 

at  values  of  and  p^  which  satisfy  the  transcendental  equation: 

& lX0^(t  ~ Q (lz± ^a) . 

P,  *>«/&/,  ’ 


2.49 


where,  fron  equations  2.48, 

„»*.  A j /■  , - % 

,«i  = # - a.(i-‘.Li)  + U 
<h*  Cl 

The  oidiho^onality  properties  of  the  F-f unctions  can  be  calculated  directly: 

H*  A 


Y'-^- 
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and  the  comolciier.tr.ry  solutions  car;  bo  chosen  as  follows: 

li?  eight  ad  orthogonality  conditions  of  the  type  found  in  the  familiar  vibrating 
iiQrcbrane  problem  uro  also  cxciuccd  here,  eince  the  pressure  field  is  not 
everywhere  twice  differentiable. 


vhare  the  final  result  contains  the  substituted  equation  2.49.  It  is  clear 
that  tho  functions  are  not  orthogonal  unless  the  corpiex  densities  are  iden- 
tical. For  the  lossless  cyst on.  the  nodes  are  orthogonal  only  in  the  absence 
of  a density  discontinuity.  This  situation  arises  from  tho  peculiar  behavior 
of  the  fiolde  at  a medium  discontinuity  as  discussed  In  Chapter  II.  The 
method  of  orthogonal  modes  is  therefore  not  generally  applicable  for  the 

two-layered  meciin,  and  we  turn  our  attention  to  the  solution  by  the  Fourier- 

1 

Bessel  transformation.-1- 

The  point  scarce  tern  in  equation  2.45a  is  again  resolved  into  a 
continuum  of  plane  sources  and  the  solution  of  the  pressure  wave  equations 
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The  pressure  fields  are  expressed  as: 


2.55b 


^ * fa*  A 2-55a 

"fo.  ""  'ffe.L  2.52b 

where  A>  ;md  A?  arc  the  particular  and  comolenentary  components  of  the 
VtP  0<G  ‘ 

total  fie  .Id  in  the  upper  medium.  The  field  in  the  lower  medium  has  no  par- 


f ’ r*n1or»  1 1 vf -4  < 
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The  particular  solution  can  be  chosen  to  bo  the  same  an  that 


derived  previously: 


fa  = r ; otifJ 


provi deu 
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The  field  of  the  point  source  can  then  be  found  in  the  sane  manner  as  before, 
and  expressed  in  a contracted  form: 

ft  -Jty  2l(kr)kPk  sinpf  2.(co$feJ)jk'  -hpjt'nft  Jsmfi  o±2±<{ 

ft  t&P  TJMMJK  *inpj[us(p,V 
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If  the  medium  were  completely  homogeneous,  the  above  ca^resentations  would 
reduce  to  the  integrals  found  for  the  case  of  the  single  homogeneous  layer, 
equations  2.55. 

The  remainder  of  this  derivation  will  be  concerned  only  with  the 
completely  lossless  medium.. 

The  above  integrals  can  be  carried  into  the  complex  Y-plane  in  the 
sam3  m<*nnor  aa  the  case  of  the  single  homogeneous  layer.  The  major  difference 
ic  ths  fcvj^v^.l.uCwnc^c  of*  ths  ini.s£i*3JK/3  sines  ihsi*s  ox  1st  four*  possibilities5 
of  sign  combinations  for  the  ancl  parameters.  We  shall  choose  the 
positive  signs  for  tho  roots  for  the  saie  reasons  discussed  before:  thus: 


> P>  * : *■**  J P = 
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To  insure  that  the  integrands  along  the  infinite  arcs  vanish  for  all  ranges 
and  depths,  excluding  the  radial  origin,  the  branch  cuts  ere  chosen  so  that 
th.o  imaginary  parts  of  both  and  pg  are  negative  over  the  entire  right 
ha3f  <£-plane.  The  cuts  are  shown  in  Figure  3. 


Figuro  8. 


The  poles  o'  tie  integrand  are  found  at  the  zeros  ol'lrj  which  is 
precisely  the  determinants!  equation  for  the  distribution  parameters  of 
3ourco-free  nodes , equation  2,49.  In  the  lossless  case,  the  transcendental 
equation  ^0  in  comparatively  easy  to  solve.  It  is  rot  necessary  to  look 
for  cc-mplex  p since  it.  is  proved  generally  lu  Part  3 that  there  cun  exist 
no  couple’/  p in  a completely  lossless  pc  dim  ir.  wh  .ten  the  fields  vanish  at 
large  ::  . The  root.s  in  the  right  half  piano  that  do  exist  fall  into  three 


40. 
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types: 


A.  A finite  number  with  Yj  real,  ‘V-g^  “ ^ tV®^»  ^1  rea^*  ^Sative  imaginary. 

B.  A finite  number  with  % real,  0 'k  < uc/cg,  Pj  and  f-2  real. 

C.  An  infinite  number  with  ^ negative  imaginary,  P-,  and  real. 


The  root3  are  demarcated  in  Figure  8 at  typical  locations. 


The  Bessel  Function  is  transformed  into  the  sum  of  the  Hankel  func- 


tions, and  the  contours  I and  II  of  Figure  9 ore  chosen.  All  of  the  line 


integrals  vanish  in  the  sane  manner  as  shown  for  the  single  layer:  only  these 


encircling  the  poles  ere  non- zero.  In  particular,  the  integrals  over  the 


branch  line  on  the  real  axis  dl3aooear  because  the  integrands  ore  even  with 


respect  to  both  p.,  and  p9.  The  final  '"on  of  the  solution  i3  then  the  simple 


sum  of  the  residues: 


*0  (j) 
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The  field  conditions  at  large  ranges  are,  of  course,  satisfied  by  the  presence 


of  the  Hankel  Functions.  The  series  converges  for  all  depths  and  ranges. 


except  for  r®0. 


The  solution  will  now  uo  discussed  in  detail. 


1.  The  power  of  the  Fourier-Bennel  integral  method  of  solution  Is  clearly 


demonstrated  in  this  problem,  since  the  non-orthogonai  character  of  the  modes 


in  the  two- layered  system  obviated  the  U3e  of  the  orthogonal  mode  technique. 


There  does  exist,  however,  a set  of  relations  among  the  modes  in  a multi- 


layered medium  which  are  analogous  to  the  ratal  orthogonality  relations 


applicable  for  the  single- layered  system.  They  are  derived  for  the  general 
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case  in  Chapter  37. 

Tho  integral  transformation  method  again  ensures  that  the  modes  form 
a complete  set  since  only  the  line  integrals  on  the  residue  contours  are 
non-vanishing. 

An  examination  of  tho  solution  reveals  that  the  modes  independently 
satisfy  all  the  boundary  conditions.  The  solution  is  also  independent  of 
ths  choice  of  branch  cut. 

2.  Ths  dependence  of  the  nodes  upon  tho  posit? on  co-ordinates  closely 
resembles  the  form  of  the  modes  in  the  single  layer  problem.  On  the  other 
hand,  the  distribution  parameters  and  have  a very  complicated  depend- 
ence upon  the  frequency  a:.<?  medium  persnot-urs,  t! trough  *hc  determ inantal 
equation  2.49.  This  tiehauior  makes  impossible  the  identification  of  the 
inodes  as  frequency  independent  entities  in  a fixed  physical  configuration. 

It  should  be  explicitly  noted  that  the  radial  dependence  of  the 
modal  field  is  identical  in  both  layers.  It  is,  therefore,  proper  that  the 
k rather  than  the  B be  termed  tho  "eigenvalues'1  of  the  wave  eouation. 

One  other  major  difference  between  the  form  of  the  modes  in  the  one- 

and  tvo-lsyerad  systems  is  the  appearance  of  an  additional  stimulation  factor, 
/_ 

> 2A/  , which, since  it  i3  a function  only  of  f?  , does  not  merit 

TTK**n  11 

further  attention.  fheso  praceeding  comments  are  also  applicable  to  the 


node  solution  for  the  next  problem. 

3.  The  similarity  of  the  modal  forms  in  the  one-  and  two-layered  cases 
allows  the  physical  interpretations  of  modes  given  on  page  11  to  be  applied 
here,  with  some  modification.  Restricting  the  discussion  to  the  completely 


lossless  medium,  we  find  that  the  freely  propagating  mode3  represent  two 
independort  forms.  The  first  form.  Type  B,  page  40.  represents  elementary 


plane  waves  with  0 < k < 


c„  which  are  propagating  at  angles  of  incidence 


less  than  critical  angle,  defined  as  & v arc  sin  c./c,,.  These  waves  travel 

^ X C. 

at  the  angle  ^ a arc  sin  k^c^/uj  in  the  upper  layer,  refract  into  the  lower 


layer  at  the  angle  &xn  " arc  3^n  k^Cg/u),  and  are  completely  reflected  from 


the  toD  and  bottom  surfaces.  Typical  ray  paths  are  given  in  Figure  9.  Tt 
ia  noted  that  the  vertical  field  distribution  i3  a purely  sinusoidal  function. 


Figure  9. 


The  second  form  of  undamped  node.  Type  A,  page  40,  represents  elementary 
plane  waves  with  Zr,<.  k < u)/c1  which  are  trapped  in  the  upper  layer  as 
shovn  in  Figure  9.  The  waves  are  propagating  at  angles  greater  than  critical 
angle , undergoing  total  reflection  at  both  the  interface  and  the  top  surface. 
Trapping  is  further  evidenced  by  the  fact  that  the  nodal  field  dependence 
upon  the  depth  co-ordinate  is  a purely  hyperbolic  function  in  the  bottom 
layer,  indicating  destructive  interference.  For  both  modal  types,  active 
power  is  transferred  only  along  the  range  co-ordinate  and  not  across  the 
horizontal  cross-section. 

Modes  whose  propagation  factors  lie  on  the  Imaginary  axis,  Type  C, 
are  purely  attenuating  radially  and  carry  no  active  energy » 
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4.  A ?yp e B node  has  cut-off  frequency  characteristics  similar  to  a mode 
in  the  single  homogeneous  layer.  Considered  a function  of  decreasing  fre- 
quency, tbs  set  of  interfering  plane  waves  travels  at  progressively  analler 
angles-  ranging  from  critical  angle  to  normal  incidence.  Below  the  cut-off 
frequency  where  0,  the  propagating  mode  disappears  and  a new  pole  on 
the  imaginary  axis  arises  as  the  nth  propagation  factor.  The  mode  is  iden- 
tified, as  the  frequency  is  changed,  only  by  its  sequence  in  the  ordered 
series  of  poles  of  the  field  integrands,  ant’  not  by  its  vertical  field  dis- 
tribution. The  Tyng  >\  node,  considered  a function  of  decreasing  frequency, 
represents  interfering  plane  waves  which  propagate  in  tho  upper  layer  at 
angles  ranging  from  grazing  incidence  to  critical  ang.lo.  'Jhen  the  frequency 
is  lowered  beyond  critical  angle,  where  k ~ cq/e,,,  the  node  disappears,  with 
no  subsequent  -appearance  of  a new  pole  on  the  imaginary  axis.  This  mode 
simply  drops  out  of  existence  beyond  the  "cut-off"  frequency?  cut-off  frequency 
for  the  T-qie  A mode  is  thus  defined  as  the  frequency  for  which  kQ=  co/c^.  The 
proof  that  the  Type  A node  cannot  exist  below  cut-off  in  the  general  case  is 
given  in  Chapter  IV. 

5.  The  ; haae  velocity-frequency  characteristic  of  the  Type  B modes  is  a 

monotonicaily  decreasing  function  fromV^*«0  at  cut-off  to  at  the  high 

frequencies.  Tho  Type  A modes  have  a similar  behavior  except  that  V**  c„ 

n c 

at  cut-off,  and  approaches  c-^  at  the  high  frequencies. 

S.  Tlie  results  can  bo  extended  qualitatively  to  '.no  general  case  of  a 
multi- lay m*ed  median  which  is  bounded  between  perfectly  reflecting  planes. 

In  all  casos,  the  field  of  an  arbitrary  source  can  be  found  in  terns  of  a 
complete  sot  of  modes*  which  are  not  orthogonal  in  the  presence  of  a density 
discontinuity.  Since  the  modes  independently  satisfy  the  source-free  wave 
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equations  and  all  the  boundary  condition?,  one  can  find  the  nodal  dependence 
on  the  position  co-ordinates  by  solving  a transcendental  equation  of  the 
type  resembling  equation  2.49.  The  amplitude  factors  must  be  found  by 
recourse  to  the  Fourior-Bes3el  integral  transformation.  In  general,  a finite 
number  of  modes  represents  waves  propagating  between  the  outer  walls,  a 
finite  number  represents  trapped  snergy  between  layers  of  large  velocity 
parameters,  and  an  infinite  number  of  strongly  attenuating  modal  fields  forma 
the  remainder  of  tho  series. 

The  Two- layered  3emi- Infinite  Medium. 

The  schematic  configuration  of  the  two- layered  semi-infinite  medium, 
which  v-'is  the  primary  concern  of  Dr.  Pekeris,  is  shown  in  Figure  10.  The 
lower  layer  represents  a liquid  bottom  which  extends  to  z~oO.  The  medium 
parameters  In  the  lower  medium  are  assumed  to  be  greater  than  the  corresponding 
parameters  In  the  upper  layer,  k point  source  is  present  in  upper. layer  at 
the  location  (0,d). 


l-H 


f>>  1 


Figure  10. 


The  major  physical  difference  between  this  problem  and  the  preceding 
ones  is  that  energy  may  leave  the  system  through  the  bottom  as  well  as  along 
the  radial  CO-  ordinate. 

The  oxact  statement  of  the  problem  is  tho  solution  of  the  pressure 


wave  equations  in  each  of  the  layers t 


?p,  = 


(l-CLt)  2.61a 


/ 


2.61b 


and  the  corresponding  velocity  fields* 


subject  to  tho  boundary  conditions* 


ield  tends  to  0 


and  r * the  fields  represent  diverging  radiation 


if  radiation  exists 


st  shown  x>hat  the  3oluuj.ons  oi 


he  source-free  wave  equation 


'9  can  easily  find  that  the  F-functions  for 


The  F-functions  are  not  orthogonal  uni 


since 


” <■  -I  I . 2.37 

The  representation  of  tho  fields  as  Fourier-Beasel  integrals  can  be 


found  by  the  previously  Established  method* 


jjh'ir.Jk 


JMfrJk.  ,1' °t 


iii'id  i a aquation  2.-33 


when  the  medium  is  completely  homogeneous,  the  last  integral  on 
right-hand  side  of  equations  2.63a  vanishes,  and  from  the  Somerfeld  tran3 


jrmation  equation  2.34,  the  acoustic  field  reduces  to  the  contributions 


in  the 


Station  2.33a  can  conveniently  be  expressed  us  the  integrals^ 

Ft  -fM  r° ( Mkdk.  s , cup,  (H-l)npS>'tf>,  (H  -ij),  d*vt  H 


In  the  abbreviated  form,  the  above  equations  read* 
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The  remainder  of  tho  analysis  deals  only  with  the  lossless  configur- 


The  Polceris  representation  of  tho  field  integrals  on  the  com  lex  2^- plane 


‘■The  same  ond  result  of  Pekeris 1 a,  with  tho  introduction  of  the  complex 
paranetera  [12] , p.  45. 


The  poles  of  the  integrand  are  those  values  of  X for  vhich  /=0. 

These  quantities  are  precisely  those  propagation  factors  of  the  modes  found 
fron  the  solution  of  tho  source-froo  wave  equation,  and  the  associated  tran^ 
seendental  aquation  2.66.  The  poles  in  the  right-half  ^-plane  fall  into  two 
categories i 

^Pekeris  [12],  Part  2. 


will  first  ba  discussed^  As  in  the  previous  problem,  there  exist  two  branch 
points  at  jj},*  pg*  0.  T ie  cuts  chosen  by  Pekeris  are  taken  parallel  to  the 
imaginary  axis  on  the  Reinann  sheet  defined  by  the  positive  signs  of  the 
roots  of  the  £ parameters,  where* 

-A  ; k'g,* " 
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A.  A finite  number  on  the  real  axle,  ; A 

r 


is  purely  x'eal,  pg  negative  imaginary. 


B.  An  infinite  number  in  the  fourth  quadrant  located  in  the  region  defined 


by  ka<  o>/cg.  Theso  polos  have  'ft  complex,  6^  and  pg  complex  uith 


positive  imaginary  parts.  It  is  particularly  noted  that  these  poles 


exist  ovon  though  the  system  is  lossless.  Dr.  Fekaris  did  not  introduce 


these  poles  in  his  formulation. 


The  poles  are  demarcated  at  typical  locations  in  Figure  10. 


Bach  of  the  integrals  of  equations  2.71  can  bo  handled  in  the  same 


manner.  Tbs  Bessel  function  13  transformed  into  the  Henkel  functions , 


referred  to  Figure  11,  the  paths  are  defined  as  follows: 


<*> 

° W ted 
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2.73a 


2.73b 


fPotes 


The  integrands  along  the  paths  CB  and  FE  vanish  as  the  radius  of  the  circular 


segments  extends  to  infinity.  The  Integra:  .d  along  the  path  HUH  vill  vanish  only 


’•hen  the  range  Is  at  sene  value  greater  than  the  depth:  this  fact  arises 


upon  examination  of  tlie  original  integrals  in  2.68  where  it  is  observed,  from 


Figure  11, that  the  imaginary  parts  of  both  the  {3  parameters  are  positive  over 
this  path.^  The  integrands  along  the  oaths  vanish  AC  and  HI  as  before.  The 


integrals  along  IIG  and  GF  cancel  because  the  integrands  are  singlo  valued 


with  respect  to  f-g  along  these  paths,  and  are  even  uith  respect  to  {3^  which 


possesses  a simple  sign  change  on  opposite  sides  of  the  cut.  The  integrals 


along  the  branch  cut  through  co/c^  do  not  cancel,  houevor,  because  the  integrands 


are  mixed  functions  of  pr>.  Tne  transformation  of  equations  2.71  then  takes  t) 

ifo  be  precise,  lm( 3]  ) and  Im(Bo)  are  positive  cn  ML,  Ifa(Bp)is  negative  and 
3ia({3^)  is  positive  on  JH. 
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The  difficulties  with  the  Pakerla  field  representation  are  observed 
to  be  these* 

1.  The  solution  is  only  applicable  for  ranges  which  are  somewhat  greater 
than  the  depth..  This  situation  is  not  only  undesirablo  from  a practical 
viewpoint,  but  it  appears  to  be  incompatible  v;ith  the  fact  that  the  solution 
for  the  hounded  systems  of  the  previous  sections  wan  valid  everywhere  Gxcept 
at  the  origin.^- 

2.  Type  D modes  contain  and  v,ith  positive  imaginary  parts.  The 
residue  of  such  a node  will  contain  the  factor  € t the  node  thus  repre- 


sents a field  which  is  exponentially 


isin/i  with  depth.  Such  peculiar 


medal  fields  are  definitely  admissable  in  the  Pekoris  formulation,  s;  ice 
the  condition  that  the  field  vanish  at  infinite  depth  is  applied  only  the 
total  field  representation.  It  should  oe  explicitly  stated  that  these  modes 
do  not  contradict  the  energy  conservation  requirements  3tated  in  Chapter  IV. 
Furthermore,  the  presence  of  a viscous  loss,  tO  thO  first  order  approxima- 
tion stated  in  Chapter  II,  will  not  cause  tho  fields  of  those  nodes  to  vanish 
at  infinite  depths. 

3.  Dr.  Pekoris  has  shown  that,  under  certain  conditions,  the  branch  line 

integral  has  the  asymptotic  form  at  large  ranges  of  a modified  dipole  field. 

At  first,  it  seers  reasonable  to  ascribe  this  physical  reality  to  the  branch 

line  integral,  since  such  a field  would  be  expected  when  the  undamped  modes 

are  all  “cut-off." However,  if  one  considers  the  Type  B modes  as  an  integral 

iFurthermore , there  is  the  ouent ion  of  the  validity  of  the  Fourier-Bessel  trans- 
formation of  the  radial  delta  function  when  the  solution  is  not  everywhere 
j>  applicable. 

In  the  bottom. 
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part  of  the  solution  (which  Pekeris  does  not),  then  it  ia  not  clear  what  a 
branch-line  integral  represents.  The  situation  is  this:  The  undamped  (Type  A) 
modes  represent  energy  trapped  in  the  upper  layer.  The  Type  B modes  represent 
tho  acoustic  energy  which  enters  the  bottom  layer,  and  comprise  an  i^iaite 
set.  All  tho  physical  expectations  are  thus  accounted  for.  What  added 
physical  feature  is  presented  by  the  non-vanishing  branch-line  integral? 

These  difficulties  are  all  removed  when  the  branch  cuts  are  taken' 
sc  that  and  pg  have  imaginary  parts  which  are  negative  over  the  entire 
i?ight  half  of  the  ^ -plane.  These  cuts  are  along  the  lines  Imagfp^J^ImagCpg]’*©, 
and  are  identical  with  the  cuts  in  Figure  3 for  the  bounded  two-layered  medium. 
Tne  only  poles  that  are  present  on  this  sheet  of  the  Reimann  surface  are  the 
Type  A poles  previously  specified.  The  Type  B pol©3  certainly  exist  but  they 
are  now  located  on  the  other  (non-physical)  sheets  of  the  Reimann  surface,  and 
consequently  do  not  appear  in  this  solution.  The  field  integrals  can  now  be 
evaluated  in  the  sane  manner  as  given  previously.  Each  of  the  integrals 
along  the  paths  of  Figure  3 vanishes,  except  those  along  the  branch  cut. 

0 < k < aye,,*  The  final  result  for  s residues  is* 


rhe  solution  converges  for  all  r and  all  z,  where  the  convergence 


of  the  branch-line  intograi  is  Justified  on  analytical  [■rounds  given 
later. 

The  solution  is  now  discussed  in  detail. 

L.  Up  are  now  able  tc-  answer  a basic  point  which  PeVeris  lest  open  for  dis- 
cussion, namely  the  reasons  for  the  existence  of  the  non- vanishing  branch- 
line integral,  and  why  it  does  not  vanish  in  the  two-layered  serai- inf inite 
nediin  which  contains  only  a velocity  discontinuity.  First  it  is  clear 
frcri  the  above  derivation  and  the  discussions  of  the  preceding  problems  that 
there  is  no  mathematical  restriction  present  which  requires  the  solution 
to  have  the  form  or  a sun  of  discrete  nodes t hence,  there  i3  no  aoriori 
reason  for  expecting  such  a solution  in  the  unbounded  m'edim,  Secondly,  the 
critical,  point  which  led  to  the  non-zero  branch-line  integral  was  the  inte- 
grand dependence  upon  the  distribution  factor  p.  In  the  semi- infinite  medium, 
it  is  recalled  that  the  integrands  were  even  with  respect  to  p^.  but  nixed 
with  respect  to  p^.  Comparing  this  fact  with  the  situation  in  the  two- 
layered  med’.tjn  vita  perfectly  reflecting  boundaries,  we  can  explain  this 
dependence  on  the  physical  basis  that  the  presence  of  a downward  traveling 
wave  alone  in  the  bottom  o'.'  the  semi- infinite  medium  appears  as  an  asymmetrical 
term  in  the  integrand.  For  a multi-layered  medium  bounded  by  perfectly  re- 
flecting planes,  there  exist  both  upward  and  downward  traveling  waves  in  each 
layer  which  are  represented  by  a symmetrical  dependence  of  the  integrand 
upon  both  the  P'3.  The  appearance  of  the  branch-line  integral  is  only  a 
consequence  of  the  fact  that  the  medium  is  unbounded  in  the  z direction,  and 
is  not  associated  with  boundary  conditions  at  the  layer  interface.  Third, 
the  necessary  presence  of  the  branch- line  integral  is  suggested  from  the 
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analogous  quantin  mechanical  problem  cf  the  one  dimensional  finite  potential 
wall.  The  solution  of  this  problem  has  the  form  of  a sum  over  a finite 


number  of  discrete  modso  (energy  lavols)  plus  an  integration  over  a continuum  i 

< 

i > 

of  "modes."  Fourth,  we  can  ascribe  a realistic  physical  interpretation  to  i 


the  branch  lina  integral,  to  be  discussed  in  a later  paragraph.  This  char- 
acterization lends  considerable  justification  for  the  validity  of  the  field 
representation  in  the  form  given  in  equations  2.75. 

2.  The  set  of  discrete  modes  is  obviously  not  complete.  As  described  later, 

these  modes  possess  cut-off  properties  similar  to  thb  Type  A modes  of  the 
previous  section.  There  exists  a frequency  below  which  all  of  the  poles  on 
the  re  if  axis  disappear,  and  no  new  poles  appear  on  the  imaginary  axis. 

An  arbitrary  field  distribution  thus  cannot  be  expressed,  solely  by  a sum  of 
discrete  modes*  the  branch  line  integral  is  present  to  account  fer  the  complete 
solution. 

It  is  thu3  possible  to  have  a solution  for  the  deter is  problem  in  which 
the  discrete  modes  arc  orthogonal  (by  eliminating  the  density  discontinuity) 
and  are  not  complete*  this  ’act  lies  in  contradistinction  to  the  solution 
for  the  case  of  the  single  homogeneous  layer. 

3.  There  exists  an  infinite  number  of  possible  representations  of  the  solu- 
tion, 3ince  the  branch  line  choice  is  essentially  arbitrary,  '-'e  know, 
however,  that  the  solution  should  be  unique.  The  difficulty  ir  resolved  by 
the  physical  interpretation  of  the  form  of  the  solution. 

It  is  recalled  that,  for  the  lossless  case,  the  fields  were  originally 
roor i seated  by  integrals  over  the  real  variable  k.  By  extending  the  path 
of  integration  into  the  complex,  plane,  we  admitted  the  possibility  of  new 
field  representations  which  are  analytically  correct  but  do  not  correspond 


to  physical  reality,  since  these  representations  contain  k as  a complex 
variable.  The  clearest  example  of  a non-physical  component  of  the  solution 
is  the  Type  D inode,  page  43,  which  is  exponentially  increasing  in  the  bottom 
layer.  The  appearance  ol  the  non-physical  forms  of  the  solution  i3  attrib- 
uted to  the  inadequacy  of  the  initial  assumptions  in  Chapter  II*  particularly 
the  neglect  of  shear  waves  in  the  fundamental  equations.  The  solution  as 
given  in  equations  2.75  is  believed  to  be  the  only  physically  meaningful 
fom,  and  in  this  sense  is  unique,  -ion-physical  solutions  for  problems  con- 
cerned with  vertically  bounded  media  are  precluded  since  the  final  forms  are 
essentially  independent  of  the  choice  of  branch  cut. 

It  should  not  be  inferred*  however,  that  the  Pekeris  solution  does 
not  havG  va?.ue  for  calculation  purposes.  Dr  Pekeric,  in  fact,  showed  an 
adequate  correspondence  between  some  numerical  integrations  of  the  original 
field  integrals,  and  tho  asymptotic  behavior  of  his  branch  line  integral. ^ 

Dr.  Ide  used  the  Type  B modes  to  illustrate  their  theoretical  correspondence 

2 

to  some  experimental  measurements  in  shallow  water.  However,  Ide  treated 

these  modes  as  though  they  were  physically  meaningful,  an  action  which  does 

not  seam  justified.  His  misinterpretation  may  have  arisen  from  an  error  in 

3 

the  choice  of  signs  in  one  of  his  introductory,  equations.  Ido’ 3 mode 
incorrectly  represents  a wave  moving  uoward  in  the  bottom  layer  with  vanishing 
amplitude  at  large  depths.  Instead  of  a wave  moving  downward  with  exponen- 
tially increasing  amplitude. 

4.  The  branch  line  integral  of  equation  2,75  is,  in  a sense,  a blend  of  the 

Type  B modes  and  tho  Pokerls  branch  line  integral  into  a physically  realistic 

^Pekeris  [12] , Figure  23. 

^Ide  [2],  App.  E. 

^Ide  [2],  App*  A,  Equation  8. 
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Ascribed  as  follows:  Ve  first  recall  that 

o taken  over  the  real  K axis,  0 </^  < -a/c.,. 


In  teres  of  the  nlane  wave  concent  of  a mode,  the  branch  line 
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sente  the  sir ’nation  of  tho  continuum  of  plane  waves  propagating  in  the  upper 
layer  ut  angles  which  are  les3  than  critical  angle,  and  undergo  multiple 
reflections  between  tho  top  surface  and  tho  interface.  Energy  ia  refracted 
into  tho  bottom  at.  each  reflection  from  the  interface.  The  individual  waves, 
represented  by  the  branch  lino  integrand  at  a particular  value  of  k,  do  not 
vanish  at  large  depths,  but  the  integrated . effect  of  these  ’laves  vanishes 
at  large  deaths  *anc  ranros)  since:'*' 

ffkUk  = £> 

wnere  fUO  is  any  reasonable  function  of  k. 

A detailed  analysis  of  the  approximate  fores  of  this  branch  line 
integral  is  beyond  the  scope  of  this  study,  Numerical  integration  nay  not 
be  too  arduous,  at  least  for  low  frequencies,  and  small  ranges  and  depths 
5.  The  trapped  nodes,  Type  A,  page  41,  have  tho  sane  general  characteristics 
as  the  Type  A.  inodes  in  the  bounded  two-layered  system.  These  nodes  represent 
sets  of  elementary  plane  waves  which  are  propagating  in  the  upper  layer  at 
discreet  angles  with  respect  to  the  vortical.  Tho  angles  fall  in  the  range 
between  between  grazing  ar.d  critical  angles.  These  waves  are,  consequently, 
totally  reflected  from  both  the  top  surface  and  the  interface.  The  node 
thus  represents  energy  trapped  in  the  upp^r  layer,  and  propagates,  without 
loss,  in  the  radial  direction.  Tho  field  of  such  a mode  is  observed  to  have 

O 

an  exponentially  decreasing  amplitude  with  depth?  and,  as  car.  be  shown  from 

the  discussioi  in  Chapter  IV,  does  not  deliver  active  energy  across  the 
lAdler  [l],  App.  4. 
xIn  the  bottom. 


CHAPTER  IV 


3r>'P  C--7.RAL  PROPERTIES  OF  THE  ACOUJPIC  MODES 


Jlltroductipn 

The  l "onaution  contained  in  this  chapter  supplements  the  discussion 
of  the  thive  examples  u:  acoustic  nodal  solutions  described  in  the  preceding 


chapter.  The  analysis  io  based  on  the  work  of  Adl<"-r  in  studies  of  electro- 

1 


magnetic  wave  propagation.  The  results  are  applicable  in  the  case  of  tb.o 
general  multi-layered  mediin  vhen  fields  have  the  form  of  eouutions  2.75  in 


which  each  co-r.oonent  va.nl shea  at  s^oc.  The  symbols  c and  v reoresent  the 

J • n — n 


pressure  ai-'  particle  velocity  of  a typical  mode.  \ branch  line  integral 

can  also  he  r-presonted  by  the  same  3jnnbol3  since  double  integrations  over  z 
and  Vf  to  te  subsequently  performed,  are  independent  operations. 

rhc  next  two  sections  of  this  part  are  devoted  to  the  derivation  of 
the  general  si t.hogonality  conditions  that  ®xist  among  the  nodes.  The  lossy 
system  is  included  insofar  an  tl.o  results  have  a simple  physical  meaning. 

The  last  section  provx fee  o physical  interpretation  of  the  propagation  factors 
for  •' lie  ^ioortta  codes  i.n  tho  multi- .layered  systems  fron  nnergy  considerations. 
dome  Useful  Iient-J ties 

Sir  es  the  modal  sum  is  a linear  c noinatior.  of  tho  solutions  of  the 
cource-frec  wave  equation,  it  is  certain  that  each  mode  is  a solution  of  the 
field  squat i or s s 


The  ri’.edi'jm  constants  are.  of  course,  those  appropriate  for  any  particular 
layer . Tho  following  Identities  c in  thereupon  be  derived  directly  from  these 
equations : 


3,1* 

5.1b 


^ ' V-  1 

y-fir,  £ = I- 

- p„*~) 

V-(p,„£  +/£%*)  = '4  L ,r  (z„-  £) 

v-(pn  Kn  rfl,  liQ  f zx„-&(-*»f>[i-<0) 

£)  ~pzi(p~.pf+p~p)  +(*£■&< 

v'fc$r  + =■  L*°t. 


3.ie 
3. Id 
3.1e 
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Power  Orthogonality  Conditions 

The  power  orthogonality  conditions  can  be  conveniently  interpreted 
as  specifications  on  the  power  outflow  due  to  nodal  interaction.  The  starting 
point  is  tho  transfer  of  the  familiar  Poynting  vector  of  electromagnetic  wave 
theory  to  the  acoustic  wave  case. 

Tho  instantaneo'is  Poynting  vector  S (r,z,t)  is  a r3ol  quantity  which 
represents  the  instantaneous  acoustic  energy  transferred/unit  area/second  at 
the  point  (r*z)  and  at  tho  timo  t,  with  a direction  indicated  by  the  space 
orientation  of  S.  From,  its  definition, 

lr  n /_ 


£ = CMp*  )l[Xe(z  e )J. 


3.2 


hh 
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he  superscript  denotes  the  complex  conjugate. 


58. 


It  is  observed  that  the  Poynting  vector  has  the  same  space  orientation  as 
the  particle  velocity  vector.  Substituting  the  complete  set  of  modes  for  the 
fields  in  the  above  equation,  vo  havei 

5 ‘[r«  e 


Ve  are  interested  in  the  total  instantaneous  energy/ second  outflow 
from  the  cylindrical  volume  V shown  in  Figure  12.  The  height  of  the  cylinder 
extends  over  the  entire  vertical  depth  of  the  medium.  The  cylinder  i3  hollow 
with  radii  and  r^  at  the  inner  and  outer  surfaces,  which  are  noted  A^  and 
Ag.  The  top  and  bottom  plana  surfaces  of  the  volume  are  noted  A^  and  A^. 
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_1_. 
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Figur 1 12. 


It  is  assumed  that  no  sources  lie  within  V.  The  inner  radius  must  be  greater 

\ 

than  zero  because  the  acoustic  field  equations  are  singular  at  the  radial  origin. 
The  total  instantaneous  power  outflow  from  V is  computed  simply  by  integrating 
the  vector  S over  the  entire  bounding  surface: 


J$'4&  = /f^  R*(p„  e‘  “'t)][%rfe  (xnR  “*)]  ‘ 4A  ■ 2-4 
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There  is  no  contribution  from  A^  and  A^  since  the  Poynting  vector  is  zero 
on  these  planes.  The  vertical  component  of  particle  velocity  is  tangential 
to  surfaces  A^  and  A^.  Equation  3.4  thus  reduces  to: 

ii  fa*  = r.  f[%  Wfte^MRek^JJSA  . w 

£ xA  i a 
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The  computation  of  the  instantaneous  real  power  transfer  is  facilitated 
when  the  modal  fields  are  expressed  as  the  following  quantities! 


Re(fne^)  = ^ 


Z 2 


3.6a 


Re(K'*eut)  - 
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3.6b 
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It  is  easily  shown  that  tha  integrands  of  equation  3.5  have  the  form: 
A.  The  n3elf-power"  term: 


gp  / _ jf’  _ it’wf'V 

-hft**r.r»e  ) 
z ~x 

3.  The  "interaction-power"  terra: 

Re  (p~  ^ t P*y"»  + P~  f/2  J . 

*.  ' 5T 

Let  us  now  consider  the  identity: 
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lossy  system 


3.13 
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In  view  of  equation  3.12  and  the  forn  of  the  power  flow  terras  of 
equations  3.7,  it  is  clear  that  the  time -varying  parts  of  the  interaction 
power  outflow  terras  vanish.  Thus  the  time-dependent  part  of  the  total  real 
power  outflow  can  be  computed  from  the  nodal  sum  as  though  the  modes  were 
propagating  independently.  For  a completely  lojsy  system,  equation  3.13  is 
validi  with  the  application  of  both  equations  3.12  and  3.13,  all  the  cross- 
product  terns  of  equations  3.7b  vanish,  and  tho  total  time  dependent  power 
outflow  i3  the  sun  cl'  trie  ;:seif- power”  toms,  equation  3.7a. 

The  tine  average  real  acoustic  power  outflow  from  the  volume  V can 
be  found  by  averaging  the  total  outflow  over  one  period  of  tho  oscillation. 
Squations  3.7  thereupon  reduce  to: 


A.  "Self-power"  terms:  fie  (fin ^m) 


3.  "Interaction-power"  terms:  dvfe  kI  +pn  kL) 


3.14a 

3.14b 


In  the  lossy  modi  in,  the  interaction  terms  give  a definite  contribution  to 
the  total  real  power  outflow.  On  the  other  hand,  these  cross-product  terms 
vanish  In  the  lossless  configuration  since  condition  5.13  is  applicable,  and 
the  time  average  outflow  car  ba  computed  from  the  sin  of  the  time  average 
"self-power"  terms,  equation  3.14a. 

The  single  homogeneous  medium  bounded  by  perfectly  reflecting  walls 
is  a special  case.  It  is  a simple  matter  to  show,  by  direct  integration  of 
the  wave  equations,  that  both  orthogonality  conditions  3.12  end  3-15  are 
valid  even  if  the  medium  is  lossy.  Thr  cross-product,  terns  in  the  power 
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outflow  sm  vanish  in  this  situation,  and  outflow  is  computed  frcm  the  simple 
sum  of  the  "self -power''  contributions. 


Enor/rv  Orthorc 


Conditions 


Vie  shall  begin  the  derivation  of  the  energy  orthogonality  conditions 
by  deriving  the  "conservation  of  energy"  theorems  for  the  acoustic  nodes. 

Let  us  consider  identity  5.1a: 

V -i,u((l'iL)tnf!£»  3.15 

This  equation  is  integrated  ovor  the  volume  in  Figure  12,  and  the  tern  on 
the  left  is  transformed  with  the  aid  of  Gauss's  Laws 


J^P^trtdA  fit  3.16 

where  p,  c,  L are  now  functions  of  z-.  The  left  hand  side  of  the  equation 


represents  the  instantaneous  "complex"  power  leaving  the  volume  through 
radiation  and  viscous  loss.  The  right  hand  side  represents  the  time  rate 
at  which  the  instantaneous  potential  and  binotic  stored  energies  are  decreasing 
within  the  volume.  Identity  3.1b  provides  a similar  relation  for  the  time 
independent  energies: 

/ a, t'J'S  i f f f tiji JVtf- £<tV?>.n 

&•«-  y y^  y. 

It  is  observed  that  the  time  independent  complex  power  outflow  through  radi- 
ation and  dissipation  is  the  difference  between  the  time  independent  potential 
and  kinetic  energy  storages. 
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of  power  orthogonality  condition  3.12.  We  hava  than: 
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An  alternative  relation  exists  uhen  tha  mth  mode  represents  a converging 


wave  where  the  sign  cf  v is  changed: 

na 


3.19 


Adding  and  subtracting  3.13  and  3.19  produces  tho  tine  dependent  energy 
orthogonality  conditions: 


e(i-iL),rsjv  *o 
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3.20a 

3.20b 


r rn  rr\ 

If  ue  assume  that  the  radii  of  the  volute  ir.  figure  12  are  separated  only 
by  an  infinitesimal  distance,  Taylor*  c:  thoorr.  s can  be  applied,  and  the  above 
equations  enter  their  final  form  as  encr.gr  storages/unit  radial  length,  _ . 


/fe* 
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, hn  l^n 

3.21a 

5.2ib 

The  time  independent  energy  orthogonality  conditions  can  be  developed 
in  the  same  -.tanner  vj’nen  tho  system  is  completely  lossless.  The  results  are 
thesa: 
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3.22b 


The  bounded  single  homogeneous  medium  i3  the  case  in  which  the  familiar 


orthogonality  relations  hold.  For  this  configuration,  it  can  be  shown  from 


54. 

an  integration  of  the  wave  equations  that  each  of  the  terms  of  equations 
5.21  and  5.22  is  independently  zero  even  if  the  medium  is  lossy. 

We  oan  new  understand  vhy  tho  nodes  in  the  two- layered  configurations 
aro  not  orthogonal  in  the  usual  3en3e.  From  3.21  and  3.22  it  is  clear  that, 
although  the  total  contribution  to  the  energy  storage  from  interacting  modes 
is  always  zero  (in  the  lossless  case),  it  is  not  necessary  that  the  potential 
and  binotic  interaction  energy  contributions  be  Independently  zero.  We  can 
expect,  however,  cn  the  basis  of  the  discussion  in  Chapter  III,  that  the 
usual  orthogonality  relations  will  hold  in  a multi-layered  configuration 
when  tne  medium  has  a uniform  "complex"  density. 

Physical  Interpretation  of  tho  Propagation  Factor  from  Energy  Considerations 
The  concern  here  is  only  with  the  discreet  modes  that  may  exist  as 
components  of  the  field  representation  in  e multi-layered  medium. 

Let  us  consider  tho  identity  3.1b  which  is  written  for  the  case 

-t*rO-'U£»*X  3.23 

pc 

This  equation  is  integrated  over  the  volume  V of  Figure  12,  and  the  left 
hand  term  transformed  by  Gauss's  Law* 


3.24 


where  the  integrals  over  tho  surfaces  and  vanish  in  view  of  the  field 
values  thereon.  f>,  c , L ore  again  functions  of  z.  Equation  3.24  can  also 
be  expressed  in  the  form: 
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Assure  that,  the  nth  mode?  i3  diverging.  At  large  ranges  where  the  asymptotic 
form  of  the  Hankel  function  is  valid,  the  fields  have  the  following  toms: 
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The  left  hand  term  of  3.25a.  can  thus  be  urittsn  as» 


3.27 


f ft r Kl^A  51  fy(*)£  ^ V/f. 
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3.29 


where  • 0(z)  is  a positive  real  quantity. 


Over  tha  particle  velocity  vector  ia  anti-parallel  with  respect  to  the 


normal  to  tha  surface.  Over  tid3  vector  is  parallel  to  tho  outward  normal. 
The  net  powar  outflow  is  then: 
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5.29 


A,* A* 


The  second  term  on  the  right  is  always  positive,  proving  that  a diverging 
cylindrical  uavo  does  act'ially  correspond  to  outgoing  radiation.  The  first 
integral  on  the  right  hand  side  of  tho  above  equation  is  always  negative. 

To  satisfy  equation  3.25a,  which  states  that  tho  energy  in  at  must  be 
greater  than  the  energy  out  at  Ag,  an  must  have  a negative  sign.  A similar 
analysis  of  tho  converging  wave  providos  the  genoral  result  that  eauation 


5.25a  is  satisfied  only  when  1c  and  a havo  oooosite  signs.  Thus  no  first 

n n 


and  third  quadrant  poles  will  ever  arise.  Since  tho  derivation  of  this 
theorem  stomed  directly  from  tho  fundamental  field  relations,  it  is  clear 
that  the  exclusion  of  waves  that  increase  in  tho  direction  of  propagation  is 
not  tho  result  of  an  extra  condition  imposed  on  the  solution. 
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The  attenuation  factor  can  be  given  30  interpretation  by  further 
manipulation  of  equations  3.25.  'Jith  reference  to  figure  12,  if  r„  i3  sep- 
arated fron  by  only  an  infinitesimal  distance,  Taylor* 3 theorem  provides 

this  approximate  fora  of  3.25a, applicable  at  largo  ranges* 

«o 

2^  = 3.30 

KtiA^* 

Thus  the  attenuation  factor  represents  the  constant  ratio  of  the  average 
viscous  power  loss/unit  radial  length  to  the  total  rawer  transferred  through 
the  radial  cross  section  at  a particular  range. 

An  interesting  incidental  point  is  the  ratio  of  real  tc  reactive 
power  in  a time  average  “self -power’*  tom  for  the  completely  homogeneous 
median.  From  previous  considerations,  this  tern  has  the  form: 


p V*  - 
f n n 


L+k»)  ufa-^L)] 


For  kn  >?°cn  the  ratio  is  approximately  l/*Lj  for  °<r\  the  ratio 

is  approximately  L.  '.Tien  Ic^0  a^.  the  ratio  is  slightly'  greater  than  1. 

The  remainder  of  this  section  is  concerned  only  with  the  purely  loss- 
less medium. 

'.I'e  shall  first  prove  the  following  theorem 

In  a lossless  multi-layered  medium,  a discreet  solution  of  the  source- 
free  wave  equation  which  vanishes  at  large  z cannot  have  a complex  propagation 
factor?  if,  in  particular,  the  medium  is  unbounded  in  the  z direction,  the 
propagation  factor  must  be  purely  real. 

Let  U3  assume  that  the  qth  mode  has  a complex  propagation  factor  ^ . 
The  mode  can  be  written  in  the  fora F,(z)R^(r)  where; 
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and,  as  usual: 
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The  velocity  naraneter  5 c new  considered  to  bo  a function  of  the 


depth,  but  In,  of  course',  constant  over  a given  layer.  The  complex  conjugates 


of  the  above  equations  are  these: 


0 f4Ff  =° 
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3.34 


1’e  next  Multiply  3.52  and  3.34  by  and  respectively,  and  integrate 


the  two  equations  over  the  entire  range  of  ?,.  The  integrands  are  not  def  ined 


at  the  values  of  z where  there  is  a discontinuity  in  the  medium  parameters, 


but  the  integrals  ncvortheloss  exist.  The  first  term  of  both  the  equations 


can  be  into  .grated  by  parts: 
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Since  the  field  is  zero  -at  f’i.j  , we  find,  after  subtraction  of  the 


above  equations,  that* 
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Sire?  the  distribution  naraneterr  2 cannot  bo  complex,  the  propagation  factor 
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X cannot  be  ccmnlex,  and  the  first  no*  c of  the  theorem  has  been  Droved, 
c ■ 


it  i 

\\  f f 


von  in  (' j i 


interface  Ixrtwf-en  too  bottom  two  Lay  era.  .»  single  node  is  assumed  to  have 
a purely  imaginary  propagation  factor,  in  vhicb  case  the  field  is  oscillatory 
vith  tine,  and  l:ac  ..n  exponential  decay  with  range.  The  real  power  transfer 
across  the  cylindri  wil  surfaces  1?  zero  as  shown  from  eouation  3.71.  The 


fields,  of  course,  are  zero  at  infinite  depth.  There  is  a non-zero  real 

power  contribution  'ran  enei  ~j  entering  tie  top  plane  surface  of  the  cylinder, 

•'-p  t 

since  the  mode  has  m dependence  in  the  bottom  layer.  The  result 

is  that  a finite  amount  of  energy  flows  into  the  cylinder,  but  no  energy 


leaves  or  is  dissipated  as  a viscous 


This  conclusion  contradicts  the 


conservation  of  sir  -gy  o sust ic.n  3.1b  . The  proof  of  the  theorem  is  now 
complete. 

for  a dircrete  node  which  exists  in  the  lossless  ;ysto_v  equation 
3.30  reads: 


Thus  either  a or  the  tine  average  integrated  cower  transferred  radially  must 

n J 

vanish.  If  a is  nero,  then,  from  the  above  theorem,  5:  b 0,  and  tie  node 
n n- 

represents  nr.  under  ped  cylindrical  rare,  - study  ol  the  acoustic  Poynting 
vector  for  this  c<r  a shown  that  such  a node  carries  cnay  real  tine  average 
power-  in  the  radio.,  direction,  'when  it  is  clear  that  only  reactive 

power  is  carried  :l>  the  field.  The  interpretation  of  a for  this  situation 
is  found  from  equ&-  Ion  3.231*  it  is  ratio  of  the  not  time  avera;  o energy 
store d/unit  radial  oogih  to  th  total,  reactive  power  transferred/ radian  at 
& given  range: 
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55.41 


'•■/e  now  turn  our  attention  to  the  real  propagation  factor  kR.  The 

• ^ 

sign  of  k » in  conjunction  with  the  time  factor  e , indicatos  the  direction 
n 

of  an  undamped  cylindrical  wave.  By  direct  computation  It  was  shown  that  the 
wave  direction  was  also  the  direction  of  radiated  energy.  However,  although 
■A  i3  the  velocity  of  the  propagating  wave,  the  velocity  of  energy  travel 
will  now  be  3houn  to  be  the  group  velocity,  0 u)/ £ k . From  identi ly  S.lgt 
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3.1E 


Let  us  integrate  this  equation  over  the  volume  7 in  Figure  12,  applying  the 
usual  boundary  conditions  over  A-  and  Ad: 

J (tZp'r*  ds*  ' ^f(p^  r p^-  <£)*  V 5-^ 

V e can  proceed  with  the  same  method  introduced  earlier.  The  field  quantities 
on  the  left  hand  side  of  the  above  equation  are  approximated  by  the  uso  of 
the  Hankel  Function  asymptotic  form.  If  this  terra  is  now  computed  at  ranges 
and  Tg  which  are  separated  by  an  infinitesimal  distance  then  there  results: 


Jgj 

dkh 


ftfnsjC.  ~h  PVjtaL  }<*  2 

'o[2f>c*  '2  " f 


The  cr011?  velocity  is  thus  interpreted  as  the  velocity  of  energy 


propagation  in  the  nth  mode  since  it  represents  the  ratio  of  the  total 


power  outflow  across  a cylindrical  cross  section  at  given  range  to  the 
energy  stored/unit  radial  length. 


REFERENCES 


R 


J 


E 


B.  Ad lor. 
Structures, 
Ide,  R. 
Frequencies 
R.  Lapwood 


Properties  of  Guided  Wave 3 on  In  omoyeneous  Cyl  j ndrlcai 
Rad.  Lab.  Elect.  Report  "o.  102,  H.I.T.  (1349). 

F.  Pout , a nd  L-  Fry,  The  Propagation  of  Sound  at  Low 

, riRL  Report  No.  3-2113  (1943). 

, "The  Disturbance  due  to  a Line  Source  in  a Semi- Infinite 


Elastic  Medium,"  Phil.  Trans.  Roy.  Doc.  Ser . A,  242.  63  (1349). 

H.  N.  Harsh,  Jr. , Theory  of  thg  Anomalous  Propagation  of  Acoustic  Waves 
in  the  Ocean .35JSL  -eport  No.  1L1  (1350). 

M.I.T.  R.L.  Sorias,  Volume  13,  Propagation  of  Short  Radio  ’■■/ayes. 
RcGrnv-Riliv  Y.  1351. 

P,  !•'.  ’tors®.  Vibration  o.nd  Sound*.  McGraw-Hill,  N.  Y.  (1948). 

P-  M.  Horse  and  H,  Feshbach,  ’Methods  of  Thooretlcal  Physic 3,  Tech. 
Press,  Cambridge  ^1946). 

National  Research  Council,  Play  sics  of  Sound  in  -the  Sea.-  Sum.  Tech.  Rep. 
Div.  6,  3 ( 194S ) . 

National  Research  Council,  Principles  of  Under /ater  Sound . Gum.  Tech. 
Ren..  Div .3,  7 1946). 


M.  t’e  ' .ids,  "The  Disturbance  due  to  a.  Line  Source  in  a Seal- Infinite 
.Elastic  Medium  with  a Single  Surface  Layer,”  Phil.  Trans.  Roy . Soc . , 

t.  VATZ  Ol  7 

• rtf.  tW»o/  j \ / • 

L.  Page,  Introduction  to  Thecreticol  Physics.  Van  Nostrand,  N.Y„  (1355) 
C.  L.  Peheris,  Theory'  of  Propagation  of  Explosive  Sound  in  Shallow 
Naves,"  Geo.  Soc.  of  Anar. , Her..  27  (1343). 

C.  L.  Pekeris,  -‘Ray  Theory'  vs.  Normal  Mode  Theory  in  Nave  Propagation 
Problems,"  Proc.  Gv.i.  .n  Am.  Hath. . 2.  Yi  (,.P!30). 


■rtial  f fercmtial  Equations  in  Hr 


ica,  Ac iconic 


cou-stic  av°  Propagation  Through  Inhano -enoous  ,*erila 


rinceton  Un.ive: 


7 (1351) 


^7005 


16.  .T.  A.  Stratton,  Elect  roriametlc  Theory,  KcGraw-Hili,  II.  Y.  (1941). 

17.  J.  L.  Wortol  and  II.  Ewin^,  '’Sxnloai'/e  Sounds  In  Shaila.  Water," 

Seo.  So? , of  Aner . ? Hen.  Z7,  (1948). 


Distribution  List  (cont'd) 


• Catholic  University  of  America 

Washington  L7,  D.  C. 

Attn i Professor  U.  F.  Horzfeld 

t 


Scripps  Institute  of  Oceanography 

La  Jolla,  California 

Attn:  Or.  Roger  Revslie  (1) 

Brown  University 

Department  of  Applied  Physics 

Providence  12,  Rhode  Island 

Attn:  Professor  R.  3.  Lindsay  (1) 

University  of  California 
Department  of  Physics 
Los  Angeles,  California 
Attn:  Profeasor  Rudnick  (l) 

Woods  Hole  Oceanographic  Institute 

Woods  Hole,  Massachusetts 

Attn:  Mr*  Vine  (l) 

Cotamanding  Officer 

Haval  Air  Development  Center 

Johnsville,  Pennsylvania  (l) 

University  of  Michigan 

ten  Arbor,  Michigan 

Attn:  Dr.  R. "ft.  Howe  (l) 

Lament  Geological  Observatory 

Torre  Cliffs 

Palisades,  Tew  York 

Attn:  Dr.  J.  Worzel  (1) 

Hudson  Laboiatory 
Columbia  University 
143  Palisades 

Dobbs  Ferry,  Hew  York,  N.  Y. 

Attn:  Dr.  E.  T«  Booth  (l) 

Department  cf  the  .iir  Force 

Cccinanding  Officer 

Air  Force  Cambridge  Re-search  Center 

230  Albany  -treat 

Cambridge,  I'aeeachisetts 

Attn:  Dr.  Norman  Haskell  (l) 

Director 

Defense  Research  Laboratory 

Univc,rs.l ty  cf  Taras 

Austin,  Texas  (1) 


...  S.-i*** 


